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ABSTRACT 

Context. AGB stars lose a large percentage of their mass in a dust-driven wind. This creates a circumstellar envelope, which can be 
studied through thermal dust emission and molecular emission lines. In the case of high mass-loss rates, this study is complicated by 
the high optical depths and the intricate coupling between gas and dust radiative transfer characteristics. An important aspect of the 
physics of gas-dust interactions is the strong influence of dust on the excitation of several molecules, including H 2 O. 

Aims. The dust and gas content of the envelope surrounding the high mass-loss rate OH/IR star OH 127.8-hO.O, as traced by Herschel 
observations, is studied, with a focus on the H 2 O content and the dust-to-gas ratio. We report detecting a large number of H 2 O vapor 
emission lines up to / = 9 in the Herschel data, for which we present the measured line strengths. 

Methods. The treatments of both gas and dust species are combined using two numerical radiative transfer codes. The method is 
illustrated for both low and high mass-loss-rate sources. Speciflcally, we discuss different ways of assessing the dust-to-gas ratio: 1) 
from the dust thermal emission spectrum and the CO molecular gas line strengths, 2) from the momentum transfer from dust to gas 
and the measured gas terminal velocity, and 3) from the determination of the required amount of dust to reproduce H 2 O lines for a 
given H 2 O vapor abundance. These three diagnostics probe different zones of the outflow, for the first time allowing an investigation 
of a possible radial dependence of the dust-to-gas ratio. 

Results. We modeled the infrared continuum and the CO and H 2 O emission lines in OH 127.8-hO.O simultaneously. We And a dust- 
mass-loss rate of (0.5 + 0.1) x 10"^ M© yr"^ and an H 2 O ice fraction of 16% + 2% with a crystalline-to-amorphous ratio of 0.8 + 0.2. 
The gas temperature structure is modeled with a power law, leading to a constant gas-mass-loss rate between 2 x 10"^ M© yr“^ and 
1 X 10“"^ M© yr“^ depending on the temperature profile. In addition, a change in mass-loss rate is needed to explain the / = 1-0 and 
/ = 2 - 1 CO lines formed in the outer wind, where the older mass-loss rate is estimated to be 1 x 10“^ M© yr"^ The dust-to-gas ratio 
found with method 1) is 0.01, accurate to within a factor of three; method 2) yields a lower limit of 0.0005; and method 3) results in an 
upper limit of 0.005. The H 2 O ice fraction leads to a minimum required H 2 O vapor abundance with respect to H 2 of (1.7 + 0.2) x 10""^. 
Finally, we report detecting 1612 MHz OH maser pumping channels in the far-infrared at 79.1, 98.7, and 162.9 pm. 

Conclusions. Abundance predictions for a stellar atmosphere in local thermodynamic equilibrium yield a twice higher H 2 O vapor 
abundance (~3 x 10""^), suggesting a 50 % freeze-out. This is considerably higher than current freeze-out predictions. Regarding the 
dust-to-gas ratio, methods 2) and 3) probe a deeper part of the envelope, while method 1) is sensitive to the outermost regions. The 
latter diagnostic yields a significantly higher dust-to-gas ratio than do the two other probes. We offer several potential explanations 
for this behavior: a clumpy outflow, a variable mass loss, or a continued dust growth. 
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1. Introduction 

Stars ascending the asymptotic giant branch (AGB) are cool 
and luminous, and they show pulsations with large periods and 
amplitudes. Their low effective temperature allows molecules 
and dust particles to form, with the dust playing an impor¬ 
tant r ole in driving the stellar wind these stars exhibit (iKwokI 
I1975I) . As such, AGB stars are important galactic factories of 
interstellar gas and dust, contributing significantly to the inter¬ 
stellar mass budget (IWhittetiri992t iTielensI l2005l) . More than 

Send offprint requests to: R. Lombaert, robinl@ster.kuleuven.be 

* Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


70 molecular sp ecies have thus far been detected in AGB stars 
(l01ofssonll2008l) . Of these, carbon monoxide (CO) is one of the 
most abundant circumstellar molecules after molecular hydro¬ 
gen (H 2 ), locking up either all carbon atoms or all oxygen atoms, 
whichever is the least abundant. When carbon is more abun¬ 
dant (i.e. the carbon-to-oxygen number abundance ratio nc/no 
> 1; defining C-type stars), the molecules and dust species 
will typically be carbon-based. When oxygen is more abun¬ 
dant (i.e. nc/no < 1; M-type), the circumstellar envelope (CSE) 
will consist niainly of oxygen-based molecule s and dust species 
(iRussellin 934 iGilmanlll 9691:iBeck et alJfT^ . 

As the star ascends the AGB, the mass loss increases grad¬ 
ually, eventually leading to the final phase, which is suggested 
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to be a superwind (lRenzini|[l98ll) . If the AGB star has not 
yet transitioned into a C-type star when it reaches the super¬ 
wind phase, it is generally known as an OH/IR star, a name that 
stems from the presence of strong hydroxyl (OH) masers and in¬ 
frared (IR) dust emission. For OH/IR stars, the comparison of 
mass-loss rates determined from the emission of low-excitation 
CO rotational transitions and those determined from the IR con¬ 
tinuum emission appear to indicate surprisingly high dust-to¬ 
gas ratios > 0.01 dHeske et al.lll99Ql: lJusttanont & Tielenll 19921: 
iDelfosse et all 1 199^ . As IR dust emission originates in re¬ 

gions closer to the stellar surface than low-excitation CO emis¬ 
sion, therefore tracing a more recent history of the mass-loss 
behavior, these high dust-to-gas ratios may be spurious and in 
reality be a manifestation of the recent onset o f a superwind 
(lJusttanont & Tielen3ll992l:lbelfosse et al.lll997l) . 

Water (H 2 O) vapor has been detected in CSEs of all chemical 
types, albeit with significantly higher abund ances with respect 
to H 2 in M-type AGB stars (^Hoo/^h. ~ 1 0“"^; ICherchnefl2006l: 
iMaercker et al.ll200^ iDecin et sd]l2010bh . In these stars, H 2 O 
vapor plays an important role in the energy balance because 
it is one of the dominant coolants in the innermost regions 
of the envelope thanks to its large number of far-IR transi¬ 
tions (iTruong-Bach et all Il999l) . It is, however, difficult to 
determine H 2 O vapor abundances accurately from H 2 O vapor 
emission, owing to, e.g., a complex ro-vibrational molecular 
structure, multiple excitation mechanisms, and s aturation effects 
(IMaercker et al.|[2QQM 120091:IPecin et ^l2QlQbl) . 

Hitherto, a lack of H 2 O observations has been hamper¬ 
ing a detailed analysis of the H 2 O excitation and abun¬ 
dance. Some H 2 O masers and vibra tionally excited H 2 O lines 
have been detected fro m the ground (iMenten & Melnicklll989l: 
iMenten et all 120061: see IMaercker et all l2008l for a summary). 
A detailed survey of multiple H 2 O vapor emission lines, how¬ 
ever, requires observations made from space. Until recently, 
only a few space missions have detected circumstellar H 2 O 
emission in the far -IR. The Infrared Space Observatory (ISO, 
iKessler et al.ll 19961) found a rich H 2 O spectrum for multiple ob¬ 
jects, though the spectral resolution was too low to detect any¬ 
thing but the strongest emission lines (|Truong-Bach et alll 19991: 
iBarlow et~n]ll996tlNeufeld et al.lll996l) . 

The recently launched Herschel Space Observatory 
(iPilbratt et al.ll201Qh . allows for a breakthrough in the study of 
H 2 O in AGB sources. OH 127.8-1-0.0 is the first high mass-loss 
OH/IR star observed with the Photodetectin^ Array Camera and 
Spectrometer (BAGS, iPoglitsch et al.l 1201^ onboard Herschel. 
High-J CO emission has also been detected in observations 
made by the Heterodyn e Instrument for the Far-Infrared (HIFI, 
Ide Graauw et al.llMl^ . We aim for a comprehensive study of 
the physics of H 2 O in OH 127.8-r0.0 by introducing a combined 
modeling of the gaseous and the solid state components of 
the outflow. We determine the gas-mass-loss rate, the radial 
abundance profile of H 2 O vapor, the location of H 2 0-ice 
formation, and the H 20 -ice characteristics, i.e. the ratio of 
amorphous to crystalline ice particles. We also address the 
dust-to-gas ratio using three different diagnostics. The first uses 
the thermal IR continuum of the dust, the second establishes the 
amount of dust needed to accelerate the outflow to the observed 
terminal gas velocity, and the third is based on the impact of 
dust emission on the strength of H 2 O lines for a given H 2 O 
vapor abundance. These three diagnostics probe different zones 
of the circumstellar envelope, for the first time allowing an 
investigation of a possible radial dependence of the dust-to-gas 
ratio. 


Table 1. Overview of some stellar and circumstellar parameters of OH 
127.8+0.0. The distance to the star is denoted as d^, the stellar lumi¬ 
nosity as L^, the CO abundance as wco/^H 2 ^ the pulsational period as P, 
the stellar velocity with respect to the local standard of rest as ulsr, the 
stochastic velocity in the outflow as t^stoch, and the gas terminal velocity 
as Voo^g • 


Input parameters 

J* 2.1 kpc 

P 1573 days 

L* l.OxlO^Lo 

vlsr -55.0 kms“^ 

WC 0 /WH 2 2.0 X 10“^* 

Voo,g 12.5 kms“^ 

«stoch 1-5 kms“' 



2. Target selection and data reduction 

2.1. The OH/IR star OH 127.8-hO.O 

OH 127.8+0.0, also known as V669 Gas, is a high mass- 
loss-rate AGB star with a relatively simple geometry. VLA 
maser maps of this object sho w an almost spherical struc¬ 
ture (iBowers & JohnstonI Il990l) . The maps hint at pos¬ 
sible dumpiness in the gaseous component of the CSE. 
Estimates for the distance to this source vary from 1.8 
kpc to 7 kpc, corresponding to a luminosity range from 


6 X 10^ L© to 2.6 X 10^ L© (Herman 

& Habind 

1985; 

Ensels et al.lll986l:lBowers & Johnstonll 19901: 

Heske et al.l 

1990; 

van Laneevelde et al.l I 1990 I: iKemper et al.l l2002l). We i 

'ollow 

:count 

Suh & KimI (l2002b. who take the pulsational phase into ac 


while modeling the spectral energy distribution (SED). They 
find a luminosity of L^,max = 3.6 x 10^ L© at light maximum, 
f-*,min = 1.0 X 10"^ L© at light minimum, and an average lu¬ 
minosity of T*,avg = 2.7 X 10"^ L©. The last agrees with th e 
period-luminosity relations derived by IWhitelock et al.l (Il99ll) . 
taking the pulsatio nal period equal to P = 1537 + 17.7 days 
(ISuh & Kiml2002l). Since the IR ISO Short Wavelength Spec¬ 
trometer (^SWS: Ide Graauw et al.lll99^ data (observed in Jan¬ 
uary 1998), as well as the PACS data (January 2010) were taken 
at light minimum, we take L* = 1.0 x 10^ L©. This value cor¬ 
responds to a distance of J* = 2. 1 kpc. We assume a CO abun¬ 
dance of nco/nu 2 = 2.0 x lO""^ (iDecin et al.ll2010ah . The gas 
terminal velocity Voo,g = 12.5 kms“^ is determined well by the 
width of the low-excitation transitions of CO (see Pig. [B, and 
is used as the primary constraint on the gas velocity field. The 
velocity of the system wi th respect to the loca l standard of rest 
is t’LSR = -55.0 kms“^ (iDe Beck et al.ll201Ql) . The stochastic 
velocity of the gas in the wind is taken to be fstoch = 1.5 kms“^ 
(ISkinner et al.lll99^ . The stellar and circumstellar parameters 
for OH 127.8+0.0 are summarized in Table[T] 


The CSE has been modeled by several authors who report 
high gas-mass-lo ss rates of Ma ~ 10~^ - 10""^ M © yr~| 


(iNetzer & Knaod 119871; Bowers & JohnstonI 

_ 

1990; 

Justtanont & TielensI119921: ILoud et al.lll993: Suh & KimI 

2002 ; 

De Beck et al.ll2010b. Owine to the hi eh mass-loss rate, the den- 


sity in the CSE is high enough for H 2 O ice t o freeze out, shown 
by a strong absorption band around 3.1 um dOmont et alll 19901: 

lJusttanont & Tielenill992l: ISvlvester et al.lll99^ 
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2.2. Observations and data reduction 

2.2.1. PACS 

We combined three PACS observations of OH 127.8+0.0 with 
six Herschel observation identifiers (obsids, 1342189956 up to 
1342189961) taken in January 2010. The first observation was 
performed with the standard Astronomical Observing Template 
(AOT) for SED. The two others were originally obtained as part 
of the AOT fine-tuning campaign. The corresponding observ¬ 
ing modes are identical to the standard one, except that alter¬ 
native chopping patterns were used. All observations were re¬ 
duced with the appropriate interactive pipeline in HIPE 8.0.1, 
with calibration set 32. The absolute fiux calibration is based on 
the calibration block (i.e. the initial part of the observation, per¬ 
formed on internal calibration sources) and is accurate to about 
20%. We have rebinned the data with an oversampling factor of 
2, i.e. a Nyquist sampling with respect to the native instrumental 
resolution. Consistency checks between the pipeline products of 
the observations made with the three chopping patterns show ex¬ 
cellent agreement, well within the calibration uncertainty. Since 
OH 127.8+0.0 is a point source, the spectrum is extracted from 
the central spaxel and then point-source-corrected in all bands. 
We list the integrated line strengths of detected emission lines 
in Table IaTT] included in the appendix. The line strengths were 
measured by fitting a Gaussian on top of a continuum to the 
lines. The reported uncertainties include the fitting uncertainty 
and the absolute fiux calibration uncertainty of 20%. Measured 
line strengths are fiagged as line blends if they fulfill at least one 
of two criteria: 1) the full width at half maximum (EWHM) of 
the fitted Gaussian is larger than the EWHM of the PACS spec¬ 
tral resolution by at least 30%, 2) multiple H 2 O transitions have a 
central wavelength within the EWHM of the fitted central wave¬ 
length of the emission line. Other molecules are not taken into 
account. We caution the reader that the reported line strengths 
not fiagged as line blends may still be affected by emission from 
other molecules or H 2 O transitions not included in our modeling. 

2.2.2. HIFI 

OH 127.8+0.0 was observed with the HIEI instrument in the 
HIEI Single Point AOT with dual-beam switching. The observed 
rotationally excited lines in the vibrational groundstate include 
the / = 5 - 4 (obsid 1342201529, observed in July 2010) and 
7 = 9-8 (obsid 1342213357, observed in January 2011) tran¬ 
sitions. These observations were made in the framework of the 
SUCCESS Herschel Guaranteed Time program (Teyssier et al., 
in prep 9 . The data were retrieved from the Herschel Science 
Archiv^ and reduced with the standard pipeline for HIEI obser¬ 
vations in HIPE 8.1. The level 2 pipeline products were then 
reduced further by first applying baseline subtraction, followed 
by the conversion to main-beam temperature with main-beam ef¬ 
ficiencies taken from the HIEI Observers’ Manual (version 2.4, 
section 5.5.2.4), and finally by taking the mean of vertical and 
horizontal polarizations. Einally, the / = 5 - 4 line was rebinned 
to a resolution of 1.3 kms“^ and the / = 9 - 8 line to a resolu¬ 
tion of 2.2 kms“^ The absolute fiux calibration uncertainty of 
HIEI is estimated to be ~ 10% according to the HIEI Observers’ 
Manual (version 2.4, section 5.7). However, owing to the low 
signal-to-noise of ~ 4 - 5 in the observed lines, we adopt a more 
conservative calibration uncertainty of 20%. 


^ http://herschel.esac.esa.int/Science_Archive.shtml 


2.2.3. Ground-based data 

Data for several rotationally excited lines of CO in the vibra¬ 
tional groundstate were obtained with the ground-based James 
Clerk Maxwell Telescope (JCMT) and the ground-based 30m 
telescope operated by the Institut de Radioastronomie Mil- 
limetrique (IRAM). The JCMT observations include the J = 
2-1 (observed in September 2002), 7 = 3-2 (July 2000), 
7 = 4-3 (April 2000) and 7 = 6-5 (November 2002) tran- 
sitions. The first th ree JCMT transitions were published by 
iKemper et al.l (l2003l). and the 7 = 6-5 tr a nsition was pre- 
sented bv iDe Beck et al.l (l2010l) . iHeske et ^ (Il990l) published 
the IRAM observations including the 7=1-0 (June 1987) and 
7 = 2-1 (Eebruary 1988) transitions. We refer to these pub¬ 
lications for the technical details concerning the data reduction. 
In this study, the 7 = 4 - 3 transition is not taken into account. 
Considering that the line formation regions of the 7 = 3-2 
and the 7 = 4 - 3 lines largely overlap, one can expect consis¬ 
tent line-integrated fluxes for the two lines when observed with 
the same telescope. No emission in the 7 = 4 - 3 observation 
is significantly detected, while a line-integrated fiux well above 
the 3(T noise level of the JCMT observation is estimated from 
the 7 = 3 - 2, as well as from the other observations. This dis¬ 
crepancy can be caused by certain model assumptions; e.g., we 
do not consider that the CO 7 = 4 level may be depopulated by 
pumping via a molecule other than CO and therefore result in a 
significantly decreased expected 7 = 4-3 emission or by an ob¬ 
servational issue, e.g., suboptimal pointing of the telescope. The 
cause of the discrepancy is not clear, so that it is safer to exclude 
the observation from the study. 

Strong CO emission at the JCMT off-source reference posi¬ 
tion contaminates the on-source 7 = 2-1 and 7 = 3-2 JCMT 
observations after background subtraction. As shown in Eig. [T] 
the lines can be fitted with an analytical function equal to the 
sum of a s oft-parabola function representing the emission profile 
(following iDe Beck et al.ll20T^ and a Gaussian function for the 
negative fiux contribution. The Gaussian component in the fit to 
both observations is centered on ~ 50 km s“^ and has a width of 
~ 1 km s“\ which is a t ypical value for the turbu lent velocity in 
the interstellar medium (iRedfield & Linskvll2008l) . assuming the 
CO emission in the off-source observation has an interstellar ori¬ 
gin. Eor the CO 7 = 2-1 and 3-2 JCMT observ ations, we use 
an ab solute fiux calibration uncertainty of 30% (iKemper et al.l 
I 2 OO 3 I) . The CO 7 = 6-5 has a low signal-to-noise ratio and 
is therefore treated as an upper limit with an absolute fiux cali¬ 
bration uncertainty of 40%. Erom the soft-parabola component 
of the 7 = 3 - 2 observation, which both has a high signal-to- 
noise and suffers less from the off-source contamination than the 
7 = 2-1 line, we derive a gas terminal velocity Voo,g ~ 12.5 
kms“^ Eor the IRAM observ ations, we use the line profiles 
published bv iHeske et al.l (Il990h . who performed a careful back¬ 
ground subtraction to avoid an off-source CO contribution. We 
assume an absolute fiux calibration uncertainty of 20% for the 
IRAM data, taking the u ncertainty involved with the background 
subtraction into account (iHeske et al.lll99Qh . 

2.2.4. Spectral energy distribution 

The SED (see Sect. 14.11) is constructed from data obtained 
by the ISO-SWS and Lon^ Waveleng th Spectrometer (LWS; 
IClegg et ^Il996t ISwinvard et al.|[T9^ instruments, as well as 
from PACS data. The ISO-SWS data were retrieved from the 
ISloan et akl (l2003l) database. The ISO-LWS data were taken 
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Fig. 1. Ground-based JCMT observations of OH 127.8-hO.O. The 
left panel shows the CO / = 2 - 1 observation in red, whereas the 
CO / = 3 - 2 is shown in the right panel. The dashed green curve 
gives a line profile fit including a soft-parabola and a Gaussian function. 
The full blue curve indicates only the soft-parabola component, which 
represents the emission coming from the CSE of OH 127.8-hO.O. The 
Gaussian component reproduces the interstellar absorption. 


from the ISO Data Archival and rescaled to the calibrated ISO- 
SWS data. The ISO-LWS data are not background-subtracted, 
whereas the PACS data are, suggesting that more flux at long 
wavelengths is expected in the ISO-LWS data owing to OH 
127.8-fO.O’s location in the galactic plane. In addition, the PACS 
photometric data at 70 yum and 160 yum (not shown here) co¬ 
incide with the PACS spectrum. The uncertainty on the ab¬ 
solute flm^_calibration_of th^^ACS photometric data is below 
15% (iGroenewegen et all l201 ll) . Taking these considerations 
into account, the ISO and PACS data agree well. The ISO-SWS 
and PACS data were all taken at the light minimum pulsational 
phase, so we assume the san ie stellar luminosity for both data 
sets and refer to the work of ISuh & KimI (l2002h for pulsation- 
ally dependent IR continuum modeling including photometric 
data. Because OH 127.8-f0.0 lies in the galactic plane, we cor- 
rected for interstellar re ddening following the extinction law of 
IChiar & TielensI (l2006h . with an extinction correc tion factor in 
the K-band of Ak = 0.24 mag (lArenou et al.lll992h . 


3. Methodology 

To get a full, consistent understanding of the entire CSE, infor¬ 
mation from both gas and dust diagnostics should be coupled. 
Kinematical, thermodynamical, and chemical information about 
the circumstellar shell is derived from the molecular emission 
lines and the dust features by making use of two radiative trans¬ 
fer codes. The non-local thermodynamic e quilibrium (NLTE) 
line ra diative transfer code GASTRoNOoM dDecin et ffll2006l 
l2010ah calculates the velocity, temperature, and density pro¬ 
files of the gas envelope, the level populations of the individual 
molecules, and the line profiles for the different transitions of 
each molecule. T he continuum radiative transfer code MCMax 
(iMin et al.ll2QQ9l) calculates the dust temperature structure and 
the IR continuum of the envelope. These numerical codes are 
briefly described in Sects. 13.1 1 and 13.21 In Sects. 13.31 to 13.51 we 
describe how the two codes are combined with an emphasis on 
the physical connections between the gaseous and dusty compo¬ 
nents. We end this section by discussing the advantage of our 
approach in light of molecular excitation mechanisms. 


and the energy balance simultaneously (iDecin et al .1120061) . We 
assume a spherically symmetric gas density distribution. The ra¬ 
dial gas velocity profile L>g(r) depends on the momentum transfer 
via collisions between gas particles and dust grains, the latter 
being exposed to radiation pressure from the cen tral sty. Thi s 
momentum coupling is assumed to be complete (iKwo 5 I1975I) . 
such that the radiative force on the dust grains can be equated 
to the gas drag force. The population of dust grains has the as¬ 
sumed size distribution 

n^(a, r) da = A(r) a~^'^ nY{{r) da, (1) 


where is the total hydrogen number density, a the radius 
of the spherical dust grain, and A(r) an abundance scale fac¬ 
tor g iving the number o f dust particles with respect to hydro¬ 
gen (iMathis et al.lll977l) . The minimum grain size considered 
is a ^m = 0.005 urn and the maximum grain size ^max = 0.25 
yum. iHofned (l2008l) suggests that large grains are needed in an 
M-type AGB CS E to be able to driv e the stellar wind through 
photon scattering. iNorris et ^ (12012h have detected these large 
grains, with sizes up to a ~ 0.3 yum, backing up our assumption 
of a maximum grain size of ~ 0.25 yum. 

The gas kinetic temperature profile Tg(r) depends on the 
heating and cooling sources in the CSE. The heating sources 
taken into account are gas-grain collisional heating, photoelec¬ 
tric heating from dust grains, heating by cosmic rays, and heat 
exchange between dust and gas. The cooling modes include 
cooling by adiabatic expansion and the emission from rotation- 
ally excited CO and H 2 O levels and vibrationally excited H 2 lev¬ 
els. As the difference between dust and gas velocity, the drift 
velocity w(a, r) directly enters the equation for collisional gas 
heating. To calculate the contribution from the heat exchange 
between gas and du st, the dust-temper ature profile T^{r) needs 
to be known as well. iDecin et ^ (l2006h approximate this profile 
by a power law of the form 





2 /( 4 +^) 


( 2 ) 


where s ^ \ (^Qlofsson in iHabing & 01ofssonll2QQ3l) . We address 
the dust temperature profile further in Sect. l3.5.fl 


3.1.2. Radiative transfer and line profiles 

The solution of the radiative transfer equations coupled to 
the rate equat ions and the calcul ation of the line profiles are 
described by IDecin et ^ (l2006l) . In this work we adopt 
MARCS theoretical model spectra (|Decin & ErikssonI l2007l: 
ICustafsson et al.|[200^ IDecin ct al.ll2010ah to improve the esti¬ 
mate of the stellar flux, as compared to a blackbody approxima¬ 
tion. This results in more realistic absolute intensity predictions 
for the less abundant molecules with stronger dipole moments 
like H 2 O, w hich are mainly excited b y near-IR radiation from the 
central star (iKnapp & Morrisiri985l) . Eor an extensive overview 
o f the moleculy data u sed in this study, we refer to the appendix 
in IDecin et al.l (l2010ah . 


3.1. Line radiative transfer with GASTRoNOoM 

3.1.1. The kinematical and thermodynamical structure 

The kinematical and thermodynamical structure of the CSE is 
calculated by solving the equations of motion of gas and dust 

^ http://iso.esac.esa.int/iso/ida/ 


3.2. Continuum radiative transfer with MCMax 

MCMax is a self-consistent radiative transfer c ode for dusty en - 
vironments based on a Monte Carlo simulation (iMin et al.l2009h . 
It predicts the dust temperature stratification and the emergent 
IR continuum of the circumstellar en velope. We use a continu¬ 
ous distribution of ellipsoids (CDE, iBohren & Huffmaiil 119831: 
iMin et al.] l2003l) to describe the optical properties of the dust 
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species. A CDE provides mass-extinction coefficients kx - or 
cross-sections per unit mass - for homogeneous particles with 
a constant volume, where the grain size (2 cde is the radius of 
a volume-equivalent sphere. The CDE particle-shape model is 
only valid in the Rayleigh limit, i.e. when A » <2 cde- For pho¬ 
tons at wavelengths A » < 2 cde , both inside and outside the grain, 
the mass-absorption coefficients /c^ are independent of parti¬ 
cle size, and the mass-scattering coefficients are negligi¬ 

ble. 

MCMax does not include a self-consistent momentum trans¬ 
fer modeling procedure, i.e. the IR continuum is calculated based 
on a predetermined dust density distribution pxir). As a standard, 
this density distribution is assumed to be smooth, following the 
equation of mass conservation M^ir) = An v^ir) pxir), with 
Mxir) = the dust-mass-loss rate, which is assumed to be 
constant, and vx(r) the dust velocity profile, which is taken to be 
constant and equal to the terminal dust velocity t>oo,d. Because 
the drift velocity is usually unknown, the dust terminal velocity 
is often equated to the gas terminal velocity Voo,g- In most cases, 
this simp lification is fo und to be inaccurate, because the drift is 
nonzero (lKwoklll9^^ . A possible improvement includes a cus¬ 
tomized density profile that takes a nonzero drift into account, as 
well as the acceleration of the dust grains derived from momen¬ 
tum transfer modeling (see Sect. 13.4.11) . In practice, the optical 
depth Ty = 1 surface in the IR lies outside the acceleration region 
for high enough dust densities, so an improved density distribu¬ 
tion in this region is not likely to affect the IR continuum of high 
mass-loss-rate stars. On the other hand, the effect on dust emis¬ 
sion features in low mass-loss-rate stars may be significant. 

3.3. The five-step modeling approach 

We solve the line radiative transfer and continuum radiative 
transfer using a five-step approach. 

1. The dust thermal IR continuum is modeled using MCMax 
to obtain an initial estimate of the dust composition, dust 
temperature, and dust-mass-loss rate. 

2. The kinematics and thermodynamics of the gas shell are cal¬ 
culated with GASTRoNOoM incorporating dust extinction 
efficiencies, grain temperatures, and the dust-mass-loss rate 
from MCMax. This provides a model for the momentum 
transfer from dust to gas, hence a dust velocity profile. 

3. Given a dust-mass-loss rate, the dust velocity profile leads to 
a new dust-density profile for which the IR continuum model 
is updated. 

4. The gas kinematical and thermodynamical structures are re¬ 
calculated with the updated dust parameters. 

5. Line radiative transfer with GASTRoNOoM is performed and 
line profiles are calculated. 

This five-step approach is repeated by changing various shell 
parameters such as the mass-loss rate and envelope sizes, until 
the CO molecular emission data are reproduced with sufficient 
accuracy. This provides a model for the thermodynamics and 
the kinematics of the envelope. The CO molecule is an excellent 
tracer for the thermodynamics of the entire gas shell because it 
is easily collisionally excited and relatively difficult to photodis- 
sociate. 

3.4. Incorporating gas diagnostics into the dust modeling 
3.4.1. Dust velocity profile 

The dust velocity profile vx(r) cannot be derived from the IR con¬ 
tinuum emission of the dust. However, the gas terminal velocity 


is determined well from the width of CO emission lines observed 
by ground-based telescopes, providing a strong constraint on the 
gas kinematical model. In conjunction with the drift velocity 
w(a, r), the gas velocity profile Vg(r) leads to vx(r). If the mo¬ 
mentum coupling between gas and dust is complete, one can 
wr ite the drift velocity at radial distance r and for grain si ze a 


as (lKwoklll975l;lTruon2-Bach et al.lll99ll;lDecin et al.ll2006h 

r 1 - 

w(a, r) = V}^(a, r) y 1 -l- x(a, r)^ - x(a, r) , with 

VYi{a,r) = . 

1 Vg(r) 

1 Qx(a)LxdA, 

x{a, r) = f 

3 

vj{r) = - 

Mg{r)c J 

1 

VY.{a,r)\ 

3kTg{ry 

pmu 

f 

2 

, and 

1/2 


Here, Qx(ci) are the dust extinction efficiencies. Lx is the 
monochromatic stellar luminosity at wavelength A, v^ir) the 
Maxwellian velocity of the gas, Tg(r) the gas kinetic tempera¬ 
ture, k Boltzmann’s constant, p the mean molecular weight of 
the gas, and mn the mass of the hydrogen atom. 

GASTRoNOoM works with grain-size dependent extinction 
efficiencies, whereas we use grain-size independent CDE mod¬ 
els for the circumstellar extinction in MCMax. As a result, the 
grain-size dependent drift velocity w{a, r) has to be converted to 
a grain-size independent average drift velocity w{r). Eor sim¬ 
plicity, we assume that the factor [ -r x(r)2 - x(r)j has a 
negligible effect. This assumption holds in the outer region of 
the CSE, where the drift velocity is much higher than the ther¬ 
mal velocity. The weighted drift velocity w(r) with respect to the 
grain-size distribution nx{a, r) from Eq.[T]can be written as 

f r) nd(a, r) da 



Assuming a grain-size distribution between lower limit a^in and 
upper limit a^^ax, this leads to 


w(r) = 


ga VK(ao, r) 
ao 


(3) 


for an arbitrary grain size a^ of a given drift velocity, with the 
weighting factor = 1.25 “ «min) ^ («max “ 

Eor GASTRoNOoM, this yields a weighting factor of ga - 0-09. 
Combining Vd(r) = w(r) Vg(r) with the equation of mass con¬ 
servation, we find a density distribution Pd(r) that can be used in 
MCMax. 


3.5. Incorporating dust diagnostics into the gas modeling 

The formation of dust species in the stellar wind has a big in¬ 
fluence on the thermal, dynamical, and radiative structure of the 
envelope; e.g., dust-gas collisions cause heating of the gas and 
drive the stellar wind, while the thermal radiation field of the dust 
is an important contributor to the excitation of several molecules, 
such as H 2 O. An accurate description of the dust characteris¬ 
tics is thus paramount in any precise prediction of the molecular 
emission. Here, we discuss the treatment of the dust tempera¬ 
ture, the inner shell radius, dust extinction efficiencies, and the 
dust-to-gas ratio. The effects of a more consistent coupling be¬ 
tween dust and gas characteristics is described in Sect. 13.61 
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3.5.1. Dust temperature and the inner shell radius 

We include an average dust-temperature profile calculated with 
MCMax in our gas modeling, instead of the power law in Eq. [2l 
This average profile is calculated assuming that the dust species 
are in thermal contact, i.e. distributing the absorbed radiation 
among all dust species such that they are at the same temper¬ 
ature at every radial point. We still use the independent dust 
temperature profiles of the different dust species - rather than the 
average profile - in the IR continuum modeling. 

The pressure-dep endent dust- c onden sation temperature is 
determined following iKama et ^ (l2009l) . setting the inner ra¬ 
dius of the dust shell. Since this inner radius indicates the 
starting point of momentum transfer from dust to gas in the CSE, 
it is assumed to be the inner radius of the GASTRoNOoM 
model as well. 


3.5.2. Dust extinction efficiencies 

iDecin et ^ (l2006h assume extinction efficiencies for spheri¬ 
cal dust particles with a dust composition typical of OH/IR 
stars, where the main component is am orphous olivine 
(Mgx,Eei_x) 2 Si 04 (lJusttanont & Tielen3ll992h . However, if one 
determines the dust composition independently by modeling the 
IR continuum, consistent extinction efficiencies can be derived. 
To convert the grain-size independent CDE mass-extinction co¬ 
efficients Kx used in MCMax to the grain-size dependent extinc¬ 
tion efficiencies Qx{o) used in GASTRoNOoM, the wavelength- 
dependent extinction coefficient is written as 

Xa = «d(a) cr^(a) = «d(a) Gi(a) ^ 

where nd(a) is the number density of the dust particles in cm"^ 
(see Eq. [T]) and crx(a) the extinction cross-section in cm^. By 
taking /cx = Xa , with px the mass density of the dust particles, 
it follows that 

4 

Qx(a) = - Kxps a, 

assuming the grains have a homogeneous grain structure. Here, 
Ps is the average specific density of a single dust grain. This con¬ 
version can be done as long as the Rayleigh assumption required 
for the CDE particle-shape mode l is v alid for every grain size a 
used in GASTRoNOoM (see Sect. l3.2l) . 


3.5.3. The dust-to-gas ratio 

The dust-to-gas ratio in AGB environments is a rather ambigu¬ 
ous quantity and is typical ly assumed to be ~ 0.005 - 0.01 
fe.g. lWhitelock et al]fl994l) . Different approaches can be used to 
estimate the dust-to-gas ratio. We assume a constant dust-to-gas 
ratio throughout the envelope in all of these definitions: 

1. Models of high-resolution observations of CO emission con¬ 
strain the gas-mass-loss rate Mg, hence the radial profile of 
the gas density Pg(r) using the equation of mass conserva¬ 
tion. The dust-mass-loss rate M^ is determined from fitting 
the thermal IR continuum of the dust. We note that the dust 
velocity field used to convert M^ into a radial dust-density 
profile Pd(r) is obtained from the GASTRoNOoM-moddmg 
and accounts for drift between dust grains and gas particles. 
The dust-to-gas ratio is then given by 


^dens 


P% 


i i \ 

\i^oo,d/ I / 


2. Given the total mass-loss rate M = Mg M^, and the 
composition and size distribution of the dust species, GAS¬ 
TRoNOoM calculates the amount of dust needed in the en¬ 
velope to accelerate the wind to its gas terminal velocity 
Voo,g by solving the momentum equation. This approach de¬ 
pends on the efficiency of the momentum coupling between 
the dust and gas components of the CSE. We assume com¬ 
plete momentum coupling, bu t we point out that this as¬ 
sumption does not always hold (iMacGregor & Stencelll 19921: 
iDecin et ^l2010ah . The empirical value of Voo,g is deter¬ 
mined from high-resolution observations of low-excitation 
emission lines, such as CO J = 1-0. The dust-to-gas ratio 
determined via this formalism will be denoted as 

3. In case of a high mass-loss rate, CO excitation is not sensi¬ 
tive to the dust emission, which allows one to constrain the 
gas kinetic temperature profile and the Mg-value by model¬ 
ing the CO emission. In contrast, a main contributor to the 
excitation of H 2 O is thermal dust emission. This allows one 
to determine the amount of dust required to reproduce the 
observed line intensities for a given H 2 O vapor abundance. 
This leads to a dust-to-gas ratio denoted as j/^h 20 , which de¬ 
pends on the adopted H 2 O vapor abundance. 

3.6. Advantages of combined dust and gas modeling: 
molecular excitation 

Calculating theoretical line profiles for molecular emission 
strongly depends on several pumping mechanisms to populate 
the different excitation levels, some of which are connected to 
the dust properties of the outflow. The most common mecha¬ 
nisms to populate the rotational levels in the vibrational ground- 
state include: 

1. Collisional excitation: Low-energy excitation is usually 
caused by collisions between a molecule and H 2 . This mech¬ 
anism becomes more important with higher densities due 
to the more frequent collisions. Eor instance, the ground- 
vibrational level of CO is easily rotationally excited (transi¬ 
tions at T > 200 pm). 

2. Excitation by the near-IR radiation field: The near-IR stellar 
continuum photons can vibrationally excite molecules. The 
vibrational de-excitation then happens to rotationally excited 
levels in lower vibrational states, with the rotational level be¬ 
ing determined by quantum-mechanical selection rules. Eor 
instance, the first vibrational state (A ~ 4.2 pm) of CO and 
the vi = 1 (T ~ 2.7 pm), y 2 = 1 (T ~ 6.3 pm), and vs = 1 
(A ~ 2.7 pm) vibrational states of H 2 O are excited this way. 
If the dust content of a CSE is high, a significant fraction 
of the stellar near-IR photons are absorbed and re-emitted at 
longer wavelengths, and cannot be used for vibrational exci¬ 
tation of molecules. 

3. Excitation by the diffuse radiation field: The diffuse field is 
mainly the result of thermal emission by dust and the in¬ 
terstellar background radiation field. These photons allow 
rotational excitation to levels that require energies that are 
too high to be accessed through collisional excitation, and 
too low to be excited by absorption from the stellar near-IR 
radiation field. Eor instance, the ground-vibrational level of 
H 2 O is rotationally excited through photons provided by the 
diffuse field (T ~ 10 - 200 pm). Increasing the dust content 
causes more pumping through this channel. 

The relative importance of these mechanisms strongly depends 
on the Einstein coefficients and on the local physical conditions 
of both the dust and gas components of the CSE. 
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. 2. The dust extinction efficiehcfesi^ivided by grain size (in cm“^) 
versus wavelength (in jj.m) used for the models shown in Figs.[3]and|4| 
At T < 25 fim the profiles are identical. From 25 //m onward, the blue 
full line and the red dashed line show a profile where the region at T > 
25 yum is replaced with a power law of the form Qext/<3 ~ assuming 
Of = 1 and Of = 2, respectively. The black full line is representative of 
a typical oxygen-rich OH/IR extinction profile as used in MCMax, for 
which the dust composition is given in Table [3] 


To show the effect of dust on line emission predictions for a 
few selected lines of CO and H 2 O in different excitation regimes 
accessible in the PACS wavelength domain, we use a standard 
input template (Table [T]) and vary one parameter at a time. We 

give an overview for high (Mg ~ 5.0 x Mg yr”*) and low 332 yhe dust opacity law 


-70 -60 -50 -40 -70 -60 

V (km s“^) V (km s“^) 

Fig. 3. Line profile predictions for the high mass-loss-rate case 
Mg = 5.0 X 10“^ Mq yr“^ The full black curve corresponds to the stan¬ 
dard model with the inner radius of the gas shell Ri g = 3 R^, the black 
extinction efficiency profile from Fig.[2|and ij/ - 0.01. In all other mod¬ 
els only a single property is modified. The dotted green curve (which 
coincides with the other curves) assumes Ri g - 10 R^, the full blue and 
dashed red curves apply the blue and red extinction efficiency profiles 
from Fig. an d the dashed-dotted magenta curve assumes if/ - 0.001 
(see Sect. l3.6l for more details). 


profile typical for an optically thick wind tracing only the termi¬ 
nal velocity. Even though dust is unimportant for the excitation 
of CO, its indirect influence through the optical depth of the in¬ 
ner region of the envelope highlights the importance of dust for¬ 
mation sequences and of the stellar effective temperature, which 
are often poorly constrained. 


(Mg ~ l.Ox 10"^ M© yr“^) mass-loss rates of the most significant 
effects including the condensation radius, the dust extinction ef¬ 
ficiency profile, and the dust-to-gas ratio. For simplicity, we as¬ 
sume a power law for the gas temperature profile corresponding 
to Model 1 in Table |4l The extinction efficiency profiles under 
consideration are shown in Fig. O We present profiles for the 
CO / = 16 - 15 transition and the H 2 O 2i,2 - lo,i and 42,3 - 4i,4 
transitions, all in the vibrational ground state. Figure [3] displays 
the high mass-loss-rate case, and Fig. |4] the low mass-loss-rate 
case. We discuss the effects below. 


3.6.1. The condensation radius 

In the high mass-loss-rate case, the condensation radius is not 
expected to have a strong influence on the theoretical line pro¬ 
files thanks to the high opacity of the envelope. Indeed, the 
full black (condensation radius /?i,g = 3 R*) and dotted green 
(Ri,g = 10 R^) models coincide in Fig. [3] and the transitions have 
a parabolic line profile typical for optically thick winds. The 
lines shown here are formed at r > 20 R* when the wind has al¬ 
ready been fully accelerated, i.e. farther from the stellar surface 
than the condensation radius used for the green model. 

In the low mass-loss-rate case, the line formation regions of 
the lines discussed here are located in the dust condensation re¬ 
gion and the acceleration zone. Increasing the condensation ra¬ 
dius in the low mass-loss rate model results in the removal of 
a relatively large amount of dust and effectively moves the ac¬ 
celeration zone outward. This manifests itself in the shape of 
the line profile. In the green model (Ri,g = 10 R*) in Fig. |4l 
the line formation regions are located where the wind is accel- 
erated. As a result, the lines exhibit a narrow G aussian profile 
(iBuiarrabal & Alcol^ll99ll: iDecin et ^l2010ah . In the black 
standard model, a narrow and a broad component are visible in 
the CO line, indicating that the line is formed both in a region 
where the wind is still being accelerated, and in a region where 
the wind has reached its terminal gas velocity. The H 2 O 2i,2- lo,i 
line, however, is only formed in the part of the wind that has just 
reached the terminal gas velocity and leans toward a parabolic 


Often, the dust extinction efficiency profile is approximated 


by 

A 


power law, 
25 jum. 


<ext 


jLamers & Cassineliil _ 

while lJusttanont & Tielens (Il992h suggest 


es pecial F 

1 Iml)' 


)ecial ly at wavelengths 
propose Of ~ 2, 

__ _ O' ~ 1 up to 1.5. 

iTielens & Allamandolal (1 19871) propose to use a ~ 2 for crys¬ 
talline grains and a ~ 1 for amorphous grains. An AGB enve¬ 
lope is usually dominated by amorp hous material (up to at least 
80 % of the dust is amorphous, e.g. lde Vries et al.ir2Q10h . How¬ 
ever, Fig. [ 2 ] shows that a = 2 is a better approximation of the 
dust extinction efficiency profile as calculated with MCMax for 
OH 127.8-rO.O. 

Comparing the three theoretical profiles for the high mass- 
loss rate case indicates the importance of the dust extinction ef¬ 
ficiency profiles. This is expected because these efficiencies de¬ 
termine the thermal emission characteristics of the grains. The 
relative change of an H 2 O line depends not only on the opac¬ 
ity law, but also on where the fine is formed in the wind and 
on the spectroscopic characteristics of the line; i.e., for different 
excitation frequencies the dust radiation field will have a differ¬ 
ent effect. It is not straightforward to predict how these changes 
will show up for given assumptions about the dust extinction effi¬ 
ciency profile. If the excitation includes channels at wavelengths 
T ~ 10 - 200 yum (i.e. excitation mechanism 3), H 2 O excitation 
is very sensitive to the properties of the dust grains in the CSE. 
At low mass-loss rates, however, the dust content is too low for 
this mechanism to contribute significantly, such that H 2 O exci¬ 
tation is controlled by the stellar radiation field in the near-IR 
(i.e. excitation mechanism 2). 


3.6.3. The dust-to-gas ratio 

At high mass-loss rates, the sensitivity of H 2 O excitation to the 
dust properties becomes very clear when comparing the low and 
high dust-to-gas ratio models in Fig.O To demonstrate this sen¬ 
sitivity, we consider first the H 2 O 42,3 “ 4i,4 fine. The excita¬ 
tion mechanism for the H 2 O 42,3 level involves first absorbing 
photons at T ~ 273 yum, where the dust radiation field is weak. 
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Table 2. Modeling results for OH 127.8+0.0, associated with Model 
2 in Table |4] gives the stellar effective temperature; and 
the dust and gas inner radii respectively; Rog the photodissociation ra¬ 
dius of ^^CO; /?o,d the dust outer radius; and Mg the dust and gas 
mass-loss rates; i;oo,d the dust terminal velocity; j/^mom, <Adens and 0 ^ h?o th e 
dust-to-gas ratios derived from three different methods (see Sect. l3.5.3l ): 
A^H 20 -ice the H 2 O ice column density; OPR the ortho-to-para H 2 O ratio; 
and WH 20 ,crit/^H 2 the critical H 2 O vapor abundance with respect to H 2 . 


and subsequently at T ~ 80 yum, where the dust radiation field 
dominates. Decreasing the dust-to-gas ratio implies that fewer 
photons are available for the channel at T ~ 80 yum, decreasing 
the population of the 42,3 level. As a result, the strength of the 
H2O 42,3 - 4 1,4 emission line is decreased significantly. Populat¬ 
ing the H 2 O 2i 2 level, on the other hand, only involves channels 
at 4 ~ 180 yum, where the dust radiation field is again weak. As 
a result, the H 2 O 2 i ,2 - lo,i line is not affected by a decrease in 
the dust-to-gas ratio. 

Both H 2 O lines are affected by a change in the dust extinc¬ 
tion efficiency profile. A profile with a different slope (cr = 1 as 
opposed to Of = 2 in this example, see Fig. O results in a re 
tively stronger dust radiation field at wavelengths 4 > 150 yum 
compared with the dust radiation field at 4 ~ 80 yum. As a 
suit, both H 2 O lines are affected because the dust radiation fi( 
becomes stronger with respect to the underlying stellar and int 
stellar background radiation field at 4 > 150 yum. CO emissior ^ 
not noticeably affected when changing the dust-to-gas ratio, ^ 
dicating that collisional excitation dominates for this molecul 

Ultimately, if collisions are not energetic enough to havi 
significant impact, it is the balance between 1) the dust, 2) 1 
stellar, and 3) the interstellar background radiation fields at 
wavelengths involved in populating a given excitation level t] 
will determine the effect of different dust properties on molecu 
line strengths. 

4. Case study: the OH/IR star OH 127.8+0.0 

We applied the combined modeling with GASTRoNOoM and 
MCMax to the OH/IR star OH 127.8+0.0. Table [2| gives the 
modeling results, which are discussed in this section. 


Modeling results 
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Fig. 5. The dust-temperature pr^Eftt^s for OH 127.8+0.0 as modeled 
with MCMax. The full colored lines indicate the specific dust species: 
cyan for amorphous silicates, red for metallic iron, blue for forsterite, 
green for enstatite, magenta for amorphous H 2 O ice, and yellow for 
crystalline H 2 O ice. Each of these profiles are cut off at the condensation 
temperature. The dashed black line gives the mean dust temperature 
profile. The full black line shows the power law from Eq. [J] with s = 1. 
The vertical dashed line indicates the inner radius of the dust shell. 


4.1. Thermal dust emission 

To model the IR continuum of OH 127.8+0.0, we followed the 
five-step approach presented in Sect. 13.31 With the assumed pa¬ 
rameters listed in Table [T] there are few parameters left to adapt 
in order to reproduce the observed IR continuum. The inner ra¬ 
dius Ri d was fixed by considering pressure-dependent dust con¬ 
densation temperatures. The stellar effective temperature 7* has 
no influence on the IR continuum of the dust due to the high op¬ 
tical depth of the wind of OH 127.8+0.0 and is constrained to 
some extent by the CO emission modeling. The dust terminal 
velocity Voo,d was derived from the momentum transfer between 
gas and dust. 

This only leaves the dust-mass-loss rate Md, the outer ra¬ 
dius of the dust shell R 04 , and the dust composition as free pa¬ 
rameters for fitting the thermal dust emission features and the 
overall shape of the IR continuum. The parameter Ro d was cho¬ 
sen such that the emergent flux at long wavelengths matches 
the PACS data well, i n agreement with the model suggested by 
iKemper et all (l2002h . ISvlvester et all (1 19991) show that the spec¬ 
tral features in the 30 to 100 yum range can be reproduced by 


a combination of amorphous sil icates, forsteri t e, ens tatite, and 
crystalline H 2 O ice. Following iKemper et al.l (l2002h . metallic 
iron was also included. The theoretical extinction coefficients 
of amorphous silicate were calculated from a combination of 
amorphous olivines with different relative magnesium and iron 
fractions, determined by modeling the d ust features in the IR 
continuum of the 0-rich AGB star Mira (Ide Vries et al.ll2QlQl) . 
The dust species and their condensation temperatures, as well as 
their mass fractions, are listed in Table [3l The dust-mass frac¬ 
tions are given in terms of mass density of the dust species with 
respect to the total dust-mass density, assuming all six modeled 
dust species have been formed. Figure [5] shows the tempera¬ 
ture profiles of each dust species. Also shown is the average 
dust temperature profile 7d,avg that is adopted as the input dust- 
temperature profile for GASTRoNOoM in our five-step approach. 
Our results for the dust composition agree well with those of 
IKemper et al.l (l2002l) . We And a higher forsterite abundance and 
slightly higher metallic iron abundance, whereas the amorphous 
silicate abundance is lower. These differences are minor. 

The mass fraction of crystalline and amorphous H 2 O ice 
is determined by fitting the 3.1 yum absorption feature in the 
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Fig. 6. The 3.1 jdm ice absorption feature. The continuum-divided 
ISO-SWS data are shown in black. The red curve gives the best fit 
model and the green curve gives the model without H 2 O ice. The dashed 
blue and dotted cyan curve give the contributions from crystalline and 
amorphous H 2 O ice, respectively. 



Fig. 7. The SED of OH 127.8-hO.O. In black the combined ISO-SWS 
and LWS data are shown; in green the PACS data are given. The dashed 
red curve is our best-fit model. The vertical dashed black line indicates 
the transition between the ISO-SWS and ISO-LWS data. 


continuum-divided ISO-SWS data, see Fig. [6l The slightly 
shifted peak position around 3.1 //m in the mass extinction coef¬ 
ficients of amorphous and crystalline ice allows one to reproduce 
the shape and strength of this absorption feature. We find a crys¬ 
talline to amorphous H 2 O ice ratio of 0.8 + 0.2 and a total relative 
mass fraction of (16 ± 2)% for H 2 O ice, which leads to a radial 
column density of A^H 20 -ice = (3.9 ± 0.5) x 10^^ cm“^. 

ISvlvester et all (Il999l) and iKemper et al.l (l2002h have mod¬ 
eled the IR continuum of OH 127.8+0.0 extensively. Using 
only crystalline H 2 O ice, they find A^H 20 -ice = 5. 5 x 10^^ cm~^ 
and ^ HoO-ice = 8.3 X 10^^ cm“^, respectively. iDijkstra et al.l 
(20061 ) have done a theoretical study of H 2 O ice formation 
( Dijkstra et all l200^ to calculate the expected H 2 O ice mass 
fractions in OH/IR stars. For a CSE with parameters similar to 
what we find for OH 127.8+0.0, they expect that only 2% of the 
total dust mass is H 2 O ice, which is a factor of 5 lower than the 
iKemper et all (l2002h results and a factor of 8 lower than our re¬ 
sults. However, they assumed an initial H 2 O vapor abundance of 
IX10“"^ in their H 2 O ice for mation models, w hich is a rather low 
estimate for an OH/IR star (ICherchnefill2006h . More H 2 O vapor 
may lead to the formation of more H 2 O ice and would be more in 
line with our results. Moreover, f ollowing their H 2 O ice forma¬ 
tion models, IDijkstra et al.l (l2006l) show that no strong H 2 O ice 
features are expected in the IR continuum at 43 yum and 62 yum 
because most of the H 2 O ice is predicted to be amorphous. Un¬ 
like this theoretical result, they point to significant fractions of 
crystalline H 2 O ice in the spectra of many sources, in agreement 
with the large crystalline fraction that we find for OH 127.8+0.0. 
They suggest several explanations for this behavior, including 
a high mass-loss rate over luminosity ratio, axisymmetric mass 
loss, and dumpiness of the wind, all of which were not taken 
into account in their ice formation models. 

For the dust composition described above we find a dust- 
mass-loss rate of = (5.0 + I.O) x 10“^ M© yr“^. This agrees 
w ell with resu lts previously obtained: = 4.0 x I0~^ M© yr~^ 

hvl^ (12004 and Mh = (7±l)xl0-^ Mq yr ^ bv iKemper et al.l 
(I 2 OO 2 I) . both assuming spherical dust grains. The use of the CDF 
particle-shape model result s in higher extin ction efficiencies rel¬ 
ative to spherical particles (iMin et al.ll2003l) . in principle imply¬ 
ing the need for less dust to fit the IR continuum of the dust. 
The choice of particle model does not significantly infiuence the 
relative mass fractions of the dust species. The resulting SED 
model, as well as the data, is shown in Eig.lTJ We lack some IR 
continuum fiux in the region 40 yum < A <10 yum in our model. 


which is a probl em that has been indic ated by previous studies o f 
OH/IR stars, e.g. IKemper et al.l(l2002l) and ide Vries et al.l(l2010h . 

4.2. Molecular emission 

We focus here on modeling the CO and H 2 O emission lines. 
Apart from these molecules, notable detections in the PACS 
spectrum concern OH emission at T ~ 79.1 yum, ~ 98.7 jum and 
~ 162.9 yum. The line strengths of these emission lines are listed 
in Table lA. Il Because the OH emission occurs in doublets, the 
li ne strengths of both co mponents have been summed. We refer 
to ISvlvester et al.l (1 1 9971) for details on OH spectroscopy. These 
detections agree with the OH rotational cascade transitions in¬ 
volved in some of the far-IR pumping mechanisi ns suggested 
as being responsible fo r the 1612 MHz OH maser (lElitzur et al.l 
Il976l:l&av k al.|[200^ . Additional OH rotational cascade tran¬ 
sitions are expected in the PACS wavelength range at T ~ 96.4 
yum and ~ 119.4 yum , but they are not detec ted. These results are 
in accordance with ISvlvester et al.l (Il997l) . who have searched 
for the 1612 MHz OH maser channels in the ISO data of the yel¬ 
low hypergiant IRC+10420. The three strongest emission lines 
were found at the same wavelengths as our OH detections, while 
the two other rotational cascade lines in the PACS wavelength 
range were significantly weaker, if detected at all. To our knowl¬ 
edge, this is the first detection of the 1612 MHz OH maser for¬ 
mation channels in the far-infrared in an AGB CSE. Owing to 
the complexity of maser formation and the spectroscopy of OH, 
however, we do not include these OH emission lines in the anal¬ 
ysis. 

4.2.1. CO emission 

We assume that the dust-to-gas momentum transfer initiates 
the stellar wind at the inner radius Ri d of the dust shell de¬ 
rived from the pressure-dependent dust condensation tempera¬ 
tures (see Sect. lOi . The outer radius Ro,g of the gas shell is 
taken as equal t o the photodissociatio n radius of CO, following 
the formalism of iMamon et al.l (Il988l) . This leaves the gas-mass- 
loss rate Mg, the stellar effective temperature T*, and the gas ki¬ 
netic temperature profile Tg(r) as free parameters to model the 
CO emission lines. 

In the five-step approach, the thermodynamics of the gas 
shell can be calculated consistently for steps 2 and 4. If the 
H 2 O vapor abundance is high (nH 2 o/^H 2 > H 2 O cool¬ 

ing becomes one of the dominant processes in the gas thermo- 
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Table 3. The dust composition of OH 127.8+0.0’s CSE. Listed are the dust species with their chemical formula, their specific density the 
condensation temperature Tcond. the mass fraction of the dust species (given as the mass density of the dust species with respect to the total dust- 
mass density Pspecies/Pd. assuming all dust species have b een formed), and the reference to the optical d ata for the opacitie s. The references for 
the op tica l constants of the dust sp ecies are as follo ws: l.lde Vries et al.l j 2010 h a nd references therein; 2 . lJager et al.l (Il998h : S. IServoin & Pirioul 
(fl^ : 4. lHenning & Stognienkol (H^ : 5. IWarredn984h : 6 . lBertie et al.Hl969h . 


Dust species 

Chemical formula 

Ps 

(g cm“3) 

Fcond 

(K) 

Pspecies/Pd 

(%) 

Reference 

Amorphous silicate 

(MgxFei_x) 2 Si 04 

3.58 

1100 

69 

1 

Enstatite 

MgSiOs 

2.80 

950 

3 

2 

Forsterite 

Mg 2 Si 04 

3.30 

950 

7 

3 

Metallic iron 

Fe 

7.87 

1150 

5 

4 

Crystalline water ice 

C-H 2 O 

1.00 

no 

7 

5 

Amorphous water ice 

a-H20 

1.00 

100 

9 

5,6 


dynamics (iDecin et alJl2006h . This introduces a significant un- 
certainty in the gas-temperature profile if the H 2 O vapor abun¬ 
dance is not well constrained. We therefore opt to parametrize 
the temperature structure. Using a grid calculation for several 
temperature structures and for a wide range of gas-mass-loss 
rates, we constrain Tg(r) empirically for OH 127.8+0.0. The 
grid probes five free parameters: the gas-mass-loss rate, rang¬ 
ing from 1.0 X 10“^ M© yr“^ to 2.0 x 10“^ M© yr“^; the stellar 
effective temperature, ranging from 2000 K to 3500 K; and the 
gas kinetic temperature profile, which is approximated by a two- 
step power law of the form Tg i(r) = T* for r < and 
= Tg \{R[) fQj. jr > We vary e\ and 62 from 0.0 
to 1.1 and the transition radius R^ from 5 to 50 R^. A power 
law with 6 = 0.5 for t he gas kinetic temp erature is expected for 
optically thin regions (iDecin et al.l120061) . but we allow for sig¬ 
nificantly steeper laws as well in view of the high optical depth 
in OH 127.8+0.0’s CSE. 

We use the spectrally resolved low-J CO transitions observed 
with JCMT and HI FI to constrain the free parameters. Following 
IDecin et ^ (l2007l) . the evaluation of the model grid is done in 
two steps. First, all models that do not agree with the absolute 
flux calibration uncertainties (Tabs on the data sets, as specified in 
Sect. O are excluded. Then, a goodness-of-fit assessment based 
on the log-likelihood function is set up to judge the shape of the 
line profile, taking statistical noise (Tstat into account. For this 
last step, a scaling factor is introduced to equalize the integrated 
intensity of the observed line profile with the integrated intensity 
of the predicted line profile. The JCMT data do not significantly 
detect the CO J = 6-5 transition. We use both the 3 (Tstat noise 
level and (Tabs to define an upper limit for the predicted intensi¬ 
ties of this line. We also compare the predicted line profiles of 
the CO / = 2 - 1 and J = 3 - 2 JCMT observations with the 
soft parabola component of the fitted line profile, rather than the 
observed line profile, in which the interstellar CO contamina¬ 
tion does not allow for a reliable determination of the integrated 
intensity and the line profile shape. 

With the exception of the CO / = 1 - 0 and / = 2 - 1 ob¬ 
servations, four models reproduce all of the available CO transi¬ 
tions, shown in Fig. [ 8 l Our estimate of the uncertainty on the 
mass-loss rates given in Table |4] amounts to a factor of three 
on the given values and is dominated by the sampling resolu¬ 
tion of the mass-loss-rate parameter in the model grid, as well 
as by the uncertainty of the CO abundance that we assume. 
These values compare well with the mass-loss-rate estimates of 
Mg ~ 5 X 10 ^ Mq yr * repor ted in the two most recent studies 
that included OH 127.8+0.0 dSiih & Kimll200l iDe Reck et all 




V (km s“^) 




Fig. 8. The spectrally resolved low-J CO observations of OH 
127.8+0.0 are shown in black. The colored curves correspond to the 
models listed in TableH) which assume a constant mass-loss rate: 1. red; 
2. blue; 3. yellow; 4. green. See Sect. l4.2.^ for further discussion of the 
validity of these CO models. 


Table 4. Values for the grid parameters of the four best fit models to the 
CO molecular emission data. Listed are the stellar effective temperature 
the powers of the 2 -step power law 6 i and 62 , the transition radius 
Rt, the gas-mass-loss rate Mg, the dust-to-gas-ratio (^dens. and the critical 
H 2 O abundance WH 20 ,crit/^H 2 - 



T 

(K) 

U 

^2 

Rt 

(R*) 

Mg 

(Moyr-i) 

^dens 

'^H20,crit/'^H2 

1 

3500 

0.2 

0.9 

5 

1.0 X 10-^ 

0.005 

8.5 X 10-5 

2 

3000 

0.2 

0.9 

5 

5.0 X 10-5 

0.01 

1.7 X 10-4 

3 

2500 

0.2 

0.9 

5 

2.0 X 10-5 

0.025 

4.0 X 10-4 

4 

2000 

0.01 

1.0 

5 

2.0 X 10-5 

0.025 

4.0 X 10-4 


l201Qh . Because the CO lines in the PACS wavelength region are 
undetected, the PACS data only provide an upper limit for the 
high-J CO emission lines. All models listed in Table |4] agree 
with this upper limit. 
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4.2.2. Validity of CO model results 

Model 1 in Table |4] requires a stellar effective temperature of 
3500 K, which is comparatively high for OH/IR stars. Ow¬ 
ing to the high optical thickness of the circumstellar shells 
in OH/IR stars, the common method of deriving stellar effec¬ 
tive temperatures based on V-K color me asurements cannot b e 
used to const r ain the effective temperature (iDe Beck et al.l201Qh . 
iLepine et"^] (1 19951) have attempted to constrain the effective 
temperatures for a large sample of OH/IR stars based on near-IR 
(K-L’) colors. They find temperatures lower than 3000 K for the 
whole sample, contrasting with the value found for our Model 
1. We choose not to exclude Model 1 because of the uncertainty 
involved in determining effective temperatures for sources with 
optically thick shells. 

All predictions in Table |4] overestimate the CO / = 2 - 1 
observations by a factor 1.5 up to 3 and the CO / = 1 - 0 line by 
a factor of 3 up to 5. Two explanations are possible: 

1. The CO / = 1 - 0 and / = 2 - 1 lines are formed in the 
outermost part of the CSE, where the contribution of the in¬ 
terstellar radiation field cannot be neglected. This radiation 
field depends strongly on the local conditions. For instance, 
if a strong UV-source is present near OH 127.8-r0.0, the pho¬ 
todissociation radi us of CO determined from the general for¬ 
malism derived bv iMamon et al.l (Il988l) would decrease. Re¬ 
ducing Ro g ~ 50 X 10^ R* to Ro,g ~ 1500 - 2000 would 
allow the model to predict the observed intensity of the CO 
/ = 1 - 0 and / = 2 - 1 lines correctly, while keeping the 
intensity of the higher-J lines the same. However, this is re¬ 
markably close to the radius o f the OH 1612 MHz maser 
shell in OH 127.8-rO.O, which iBowers & Johnston! (Il990l) 
found to be (1.38 + 0.14)". This translates to ron ~ 1000 - 
2000 Ri, at a distance of 2.1 kpc, depending on the assumed 
temperature at the stellar surface. This suggests that such a 
small outer CSE radius is unlikely for OH 127.8-hO.O. 

2. The mass loss in OH 127.8-rO.O ma y be variable, as sug - 
gested by several previous studies fe.g. lDe Beck et al.ll2010l) . 
If the mass-loss rate has been lower in the past, then the low- 
J lines might have a lower intensity compared to our pre¬ 
dictions assuming a constant mass-loss rate. To improve the 
prediction of the / = 1-0 and 7 = 2-1 CO lines, we cal¬ 
culated models with a change in mass-loss rate going from 
Mg = 1 X 10“^ M© yr“^ in the outer wind up to Mg as listed 
in Table|4|for the inner wind. The transition from high to low 
mass-loss rate occurs gradually at the radial distance Ryu of 
~ 2 500 - 4000 which translates to ~ 7.5 - 14.5 x 10^^ 
cm. iDelfosse et all (1 19971) found similar results based on the 
IRAM ^^CO and ^^CO 7 = 2-1 and 7=1-0 transitions with 
an older, low mass-loss rate of Mgj < 5 x 10“^ M© yr“^ and 
a recent, high Mg,h between 5 x 10“^ and 5 x 10“"^ M© yr“^ 
They found a transitional radius of Ryu ~ 1.8 - 5.3 x 10^^ 
cm, depending on Mg,h- Our estimate of Ryu is larger, but 
we have a stronger constraint on Ryu due to the higher-J 
CO transitions. The values we find for Ryu translate to an 
increase in the mass-loss rate in OH 127.8-r0.0 in the last 
2000 up to 4000 years, depending on Mg^h and the tempera¬ 
ture structure. This recent change in mass-loss rate is com¬ 
monly referred to as the recent onset of the superwind, which 
is often suggested for many OH/IR stars by sev eral studies 
(lJusttanont & Tielen3ll992L IDelfosse et al.ll 19971 de Vries et 
al. in prep, Justtanont et al. in prep). 

The assumption of a change in mass-loss rate to predict the low- 
J CO line strengths correctly does not affect further modeling of 


other emission lines, as long as these lines originate in a region 
within the radial distance Ryu- This is the case for the H 2 O 
vapor emission lines detected in the PACS wavelength range, so 
we use the four models listed in Table |4| in what follows. 

4.2.3. H 2 O emission 

To determine the H 2 O vapor abundance, we use j/^dens and adopt 
the gas kinetic temperature law and gas-mass-loss rate of Model 
2 in Table|4]because the mass-loss rate is closest to the estimates 
of previous studies. What follows has been done for every model 
in Table IH and even though the resulting values scale with the 
mass-loss rate, the general conclusions do not change. 

We have selected 18 mostly unblended, non-masing H 2 O 
emission lines in the PACS spectrum to fit the GAS TRoN OoM 
models. The selection of lines is indicated in Tab lelA.il We 
assume an ortho-to-para H 2 O ratio {OPR) of 3 (iDecin et al.l 
l2010bh . When using j/^dens = 0.01 derived from fitting CO emis¬ 
sion and the thermal IR continuum (see Sect. 13.5.3b for Model 2 
in Table m we find an unexpectedly low H 2 O vapor abundanc^l 
^H 2 o/^H 2 ~ 5 X 10“^, as com pared with ^Hoo/^ H o ~ 3 x 10""^ de¬ 
rived f rom chemical models (ICherchneflll2006l) . iMaercker et al.l 
(I 2 OO 8 I) also found an H 2 O vapor abundance of ~ 10“^ for the 
OH/IR source WX Psc, indicating that such a discrepancy has 
been found before in sources that have a high mass-loss rate. 

To resolve this discrepancy, we determine for a wide 
range of H 2 O vapor abundances such that our model reproduces 
the H 2 O emission spectrum of OH 127.8-1-0.0. The results for 
Model 2 in Table |4| are shown in Fig. [9] and give further clues to 
the excitation mechanism of H 2 O vapor in the high mass-loss- 
rate case. At values > 10“^, 0 ^h 2O correlates with the H 2 O vapor 
abundance. Here, pumping through excitation by the dust radi¬ 
ation field plays an important role. For lower dust-to-gas ratios, 
the dust radiation field becomes negligible for H 2 O vapor excita¬ 
tion causing the correlation between and ^h 2 o/'^H 2 to level 
off. The correlation between j/rH 20 and the H 2 O vapor abundance 
depends on the gas-mass-loss rate. For comparison, equivalent 
results for Model 1 in Tableware shown in Fig.[9l 

Figures [T0| and [TT] show the continuum-subtracted PACS 
spectrum compared to the predictions of Model 2 in Table |4] 
for ^h 2 o/^H 2 = 3 X 10“"^ and 0 ^h 2 O = 0.003. Included and in¬ 
dicated on the spectrum are all rotational transitions in the 
vibrational groundstate and all 0 -H 2 O and P-H 2 O transitions in 
the vibrational groundstate and the vi = 1 and y 2 = 1 vibra¬ 
tional states with rotational quantum number up to 7upper = 8 
in the PACS wavelength range, regardless of being detected or 
not. The 18 H 2 O transitions used in the initial fitting procedure 
are indicated as well. We calculated model spectra for the other 
temperature and density profiles in Tableland arrive at the same 
overall result as for Model 2 with some small differences in the 
relative line strengths of the lines. 

4.2.4. Validity of H 2 O model results 

A slight downward trend is present in the comparison between 
model predictions and the observations, as shown in Figs. M 
and [m with a systematic overestimation at short wavelengths 
and a systematic underestimation at longer wavelengths. This 
difference is within the 30% absolute uncertainty calibration of 
the PACS data. However, a relative trend between short and 

^ H 2 O vapor abundances are always given for ortho-H 20 alone, while 
H 2 O column densities and H 2 O ice abundanees always inelude both 
ortho- and para-H 20 . 
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Fig. 9. OH 127.8+0.0 H 2 O emission spectrum modeling results for 
the temperature law and mass-loss rate of Models 1 and 2 in Table |4] 
in red and black, respectively. j^h 20 and its uncertainty is determined 
for a wide range of (ortho + para) H 2 O vapor abundances. From the 
modeling of the IR continuum and the CO data, a value of j^dens = 0-01 is 
determined. The expected H 2 O vapor abundance from chemical models 
is 3 X 10""^ (ICherchneffll2006h . Both values are indicated by the dashed 
black lines. The dark gray area indicates the lower limit defined by 
the critical H 2 O va por ab undance derived from the H 2 O ice fraction of 
Model 2, see Sect. 14.2.51 For comparison, the light gray area indicates 
the lower limit found for Model 1. 


long wavelengths in the model-to-data comparison is unexpected 
from the absolute calibration errors. A relative uncertainty be¬ 
tween short and long wavelengths can be caused by pointing er¬ 
rors of the telescope, but this effect is likely too small to explain 
the trend that we find. This trend is present for all models in Ta¬ 
ble m although less evident for Models 3 and 4 (with the lower 
mass-loss-rate estimate of Mg ~ 2 x 10“^ M© yr“^). 

Based on this, one could opt to exclude Models 1 and 2. 
However, H 2 O is not a good tracer of the density and temperature 
structure owing to the complexity of H 2 O excitation mechanisms 
and possible maser effects. Normally, CO is a good density and 
temperature tracer, because CO is dominated by collisional ex¬ 
citation and does not mase. However, for OH 127.8+0.0, CO 
lines are optically thick and were not reliably detected in the 
PACS wavelength range. In this case, ^^CO lines would be a bet¬ 
ter tracer, but they are significantly weaker than ^^CO emission 
lines and, as such, are not detected at all in the PACS observa¬ 
tions. As long as the majority of H 2 O lines are reproduced well 
over a wide range of wavelengths in the PACS data, for which 
the signal-to-noise ratio is low especially at short wavelengths, 
we consider a model to be satisfying. Thus, we choose not to 
exclude any models based on the trend in the predictions. 

This large a set of H 2 O lines has not been modeled before in 
such detail, covering full radiative transfer modeling of the CSE 
of a high mass-loss-rate OH/IR star. The consistent prediction 
of line-integrated fluxes of H 2 O lines across a wide wavelength 
range that is well within the absolute flux calibration of the PACS 
instrument - especially in the red bands - is remarkable, consid¬ 
ering the large number of H 2 O lines and the complexity of the 
problem. 

4.2.5. The H 2 O vapor-ice connection 

An additional constraint can be placed on the estimate of the 
H 2 O vapor abundance and the associated 0 ^h 2 O- The presence of 
H 2 O ice in an OH/IR CSE provides a lower limit on the H 2 O 
vapor abundance. The con densation temperatu re of H 2 O ice is 
Tcond,ice = 110 K, following I Kama et ^ d2009l) . The condensa¬ 
tion radius associated with 110 K is Rcond,ice = 1.2x 10^^ cm. The 


Fig. 12. Schematic representation of the CO (red), the ortho-H20 
(blue), and para-H20 (green) abundanee profiles. The vertieal dashed 
blaek line indieates the H 2 O iee eondensation radius. The vertieal full 
black line indicates the location of the OH 1612 MHz maser shell, as¬ 
suming a distance of 2100 pe. The signal-to-noise of the PACS is too 
low to trace the drop in H 2 O vapor abundance (shown here for a freeze- 
out of ~ 40%) at the H 2 O iee condensation radius. 


line formation region for all unblended, nonmasing H 2 O vapor 
lines in the spectrum of OH 127.8+0.0 is mostly within this ra¬ 
dius. We can therefore define a critical H 2 O vapor abundance at 
r < Rconddce, below which there would not be enough H 2 O vapor 
to form the observed amount of H 2 O ice at Rcond,ice- Following 
our modeling of the H 2 O ice feature, the H 2 O ice column density 
at r > Rcond,ice is 8.3 X 10^^ cm“^, which leads to a critical (ortho 
+ para) H 2 O vapor abundance of ^H 20 ,crit/'^H 2 = (1 -7+ 0 . 2 )x 10 “^. 
This critical abundance depends on the gas-mass-loss rate, as the 
ice mass is compared to the equivalent molecular hydrogen mass 
in the ice shell. As such, this critical value will be different for 
the three mass-loss rates given in Table IH For comparison, the 
critical H 2 O vapor abundance is shown in Fig p for Models 1 
and 2 in Table (H The work by iDijkstra et alJ (l2003l) suggests 
that at most 20 % of the H 2 O vapor will freeze out onto dust 
grains. The actual H 2 O vapor abundance is thus expected to be 
larger than the critical H 2 O vapor abundance. Figure [12] gives a 
schematic representation of what the H 2 O vapor abundance pro¬ 
file might look like, taking H 2 O ice condensation into account 
(with a freeze-out of ~ 40%, a value that is arbitrarily chosen) 
and H 2 O vapor photodissociation in the outer envelope. 

Observational evidence of a larger actual H 2 O vapor abun¬ 
dance than the critical H 2 O vapor abundance is given by the 
presence of emission from the OH maser at 1612 MHz in a shell 
at r > Rcond,ice- The photodissociation of H 2 O into OH and H is 
one of the main OH production paths, throughout the whole en¬ 
velope, as long as i nterstellar UV radiation is available to break 
up H 2 O molecules. iNetzer & Knapp! (Il987l) have shown that the 
OH abundance reaches a maximum at the radial distance where 
the OH maser shell at 1612 MHz is observed, indicating that 
other methods of OH production closer to the star, such as shock 
chemistry, can be ignored. As a result, H 2 O needs to be present 
in the CSE at least up to the radial distance where the OH abun¬ 
dance peaks. In the case of OH 127.8+0. 0, the radius of the 
OH 1 612 MHz maser shell is (1.38+0.14)" (iBowers & Johnston! 
I 1990 I) . which tra nslates to rpH = (4.3 + 0.4) x 10^^ cm at a dis¬ 
tance of 2.1 kpc. iNetzer & Knapp! (Il987l) also give a formula for 
the expected OH 1612 MHz maser shell ra dius, which de pends 
on the assumed interstellar radiation field (lHabinglll968h . As¬ 
suming the average Habing field, we find 6.1 x 10^^ cm, whereas 
the high Habing field leads to a shell radius of 4.3 x 10^^ cm. 

Our results for the critical H 2 O vapor abund a nce ag ree well 
with those found in other studies. ICherchn^ (l2006l) derived 
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Fig. 10. The continuum-subtracted PACS spectrum of OH 127.8-hO.O is shown in red for the blue bands. The PACS band is indicated in the 
lower left corner of each spectrum. Model 2 in Table |4| with wh 2 o/^H 2 = 3 x 10""^ and j^h 20 = 0.003 is given in blue. The other parameters are 
listed in Tables [T] and [2| The colored vertical lines indicate the molecule contributing at that specific wavelength, with full black for ^^CO, yellow 
for ortho-H20, and cyan for para-H20. The dashed black-colored lines indicate the water lines used for the initial H 2 O line fitting. The forsterite 
feature at ~ 69 fim (not completely removed during continuum subtraction) and the OH rotational cascade lines ^ni/ 2 / = 1/2 n 3 / 2 / = 3/2 and 

^ni/ 2 / = 5/2 ni/ 2 / = 3/2 at ~ 79.1 jjm and ~ 98.7 jjm, respectively, (not included in our modeling) are labeled. 
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Fig. 11. The continuum-subtracted PACS spectrum of OH 127 . 8 - 1 - 0.0 is shown in red for the red bands. The PACS band is indicated in the lower 
left corner of each spectrum. Model 2 in Table | 4 ] with wh2o/^H2 = 3 x 10 “"^ and j/rH20 = 0.003 is given in blue. The other parameters are listed 
in Tables [T] and |2] The colored vertical lines indicate the molecule contributing at that specific wavelength, with full black for ^^CO, yellow for 
ortho-H20, and cyan for para-H20. The dashed black-colored lines indicate the water lines used for the initial H 2 O line fitting. The OH rotational 
cascade line ^ni/2/ = 3/2 ni/2/ = 1/2 at ~ 162.9 jjm (not included in our modeling) is labeled. At ~ 144.9 yum, another strong line appears 

both in band RIA and band RIB. A different wavelength sampling causes the line in band RIB to appear weaker, but the integrated line fluxes 
of both lines are within the absolute flux calibration uncertainty of PACS. This line remains unidentified. The CO line alone cannot explain the 
observed integrated line flux at this wavelength. 
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the expected abundances for several molecules from thermo¬ 
dynamic equilibrium and shock-induced NLTE chemistry, and 
found ^h 2 o/^H 2 ~ 3.0 x 10""^ in 0-rich AGB stars. H 2 O vapor 
abundances derived by iMaercker et al.l (l2008l) for most sources 
in their sample lie between /^H 2 = 2.0 x 10“^ and 1.5 x 10“^. 
They do find a remarkably low H 2 O vapor abundance of ~ 10“^ 
for the OH/IR star WX Psc and offer two explanations: 1) H 2 O 
ice formation depletes H 2 O in gaseous form, and 2 ) H 2 O lines 
may be forme d in a region of a mor e recent, lower mass-loss 
rate. However. IMaercker et al.l (l2008l) use j/^dens for their molec¬ 
ular emission modeling. Our modeling has indicated that j/^dens 
leads to too low an H 2 O abundance in OH 127.8-rO.O when com¬ 
pared to the H 2 O ice content. Given that OH 127.8-rO.O and WX 
Psc are similar, their low value for the H 2 O vapor abundance in 
WX Psc could also be the result of the use of </^dens as an estimate 
of the dust-to-gas ratio. 

4.3. Discussion: The dust-to-gas ratio 

The dust -to-gas ratio typical of A GB circumstellar environments 
is 0.005 (IWhitelock et al.|[T9^ . We derive different values de¬ 
pending on the method used (for Model 2 in Table ID): 

1. We find j/^dens = 0.01, accurate to within a factor of three, 
from IR continuum and CO molecular emission modeling. 

2. The momentum transfer equation leads to a dust-to-gas ratio 
<Amom = 0.0005, while assuming complete momentum cou¬ 
pling between dust and gas, i.e. dust grains of every grain 
size are coupled to the gas. For circ umstellar environments 
typica l of stars like OH 127.8-r0.0, iMacGregor & Sten^ 
(Il992h find that silicate grains with initial size smaller than 
~ 0.05 Jim decouple from the gas near the condensation ra¬ 
dius. If the coupling is not complete, a higher dust content 
is required to arrive at the same kinematical structure of the 
envelope, implying that j/^mom > 0.0005. 

3. The critical H 2 O vapor abundance provides a strong con¬ 
straint on the expected initial H 2 O vapor abundance in OH 
127.8-r0.0. Applying our value of ^H 20 ,crit/^H 2 > 1.7 x 10“"^ 
to the grid calculation shown in Fig.[9l we find an upper limit 
for the associated dust-to-gas ratio i//u 20 < 0.005. This upper 
limit takes the uncertainty shown in Fig. [9] into account and 
assumes the unlikely case of 100 % freeze-out of H 2 O vapor 
into H 2 O ice. For reference, assuming a freeze-out of 20%, 
we arrive at 0 ^H 2 O-fo ~ 0.0015, accurate to within a factor of 
two. 

The results obtained for the dust-to-gas ratio appear incom¬ 
patible. However, each method traces a different part of the en¬ 
velope (see Fig. [13. 

1 . <Adens is based on modeling the thermal dust emission, which 
traces the dust content of the envelope out to a radius of 
~ 5000 R*, and the CO / = 9 - 8 down to / = 3 - 2 emission 
lines. These lines are formed in the outer regions of the CSE 
at 100 R* < r < 4000 R*, for Model 2 in Table (H Assum¬ 
ing the dust-mass-loss-rate remains constant throughout the 
whole CSE, j/^dens is therefore sensitive only to the outermost 
region of the envelope. 

2 . <Amom is determined from the momentum transfer equation 
and therefore traces the acceleration zone, which in our 
model is located at r < 50 R*. 

3. <Ah 20 traces the outflow at 20 R* < r < 800 R* where all 
of the H 2 O emission lines used to determine the H 2 O abun¬ 
dance are formed. 



Fig. 13. The results of determining the dust-to-gas ratio using the 
three different methods described in Sect. 13.5.^ are shown for Model 2 
in Table |4] The horizontal bar indicates the part of the envelope traced 
by the method. The vertical bar indicates the uncertainty on the indi¬ 
cated value, i^rnom is a lower limit, whereas j^h 20 is an upper limit. For 
reference, the dashed red line indicates j/^H 20 -fo assuming 20% freeze- 
out of H 2 O vapor into H 2 O ice. The relative differences between the 
three values of the dust-to-gas ratio for the other models in Tableware 
similar, but scale upward or d ownward uniformly depending on the gas- 
mass-loss rate. See Sect. 14.3I for a more detailed description. 


As shown in Fig. [m our findings tentatively point to the pres¬ 
ence of a gradient in the dust-to-gas ratio with radial distance. 
The results shown here are for Model 2 in Table IH but the same 
relative differences between the dust-to-gas ratio estimates are 
seen for the other models in Table [U There could be several 
potential explanations for such behavior. 

First, we assume a constant mass-loss rate for both the gas 
and dust components. If OH 127.8-fO.O’s mass-loss history is not 
constant, a recent increase in the gas-mass-loss rate can explain 
the gradient in the dust-to-gas ratio only when the dust-mass- 
loss rate has not increased by the same factor as well. This is 
possible only if the dust forms less efficiently for an increased 
gas density. There is no immediate evidence that suggests such 
behavior for higher mass-loss rates, so this scenario appears to 
be unlikely. 

Second, 84% of the dust mass is formed in the innermost 
region of the envelope, at a few stellar radii, and the dust-mass- 
loss rate is assumed to be constant. If dust formation extends 
beyond the vicinity of the dust condensation radius, this could 
explain the gradient in the dust-to-gas ratio. H 2 O ice formation 
is possible at a radius of ~ 1000 in the case of Model 2 in 
Table m owing to the high H 2 O vapor abundance. However, the 
amount of H 2 O ice formed is not enough to explain the radial 
increase in dust-to-gas ratio. Formation of other dust species 
(such as silicates) at large distances from the star is unlikely due 
to the lower densities of the precursor molecular species when 
compared to H 2 O vapor. 

Third, we do not take clumping into account in the models. 
If clumps are present in the envelope, the ones close to the stellar 
surface are likely to be much more optically thick than those in 
the outer envelope. As a result, we trace the real amount of gas 
and dust in the outer envelope, whereas we trace a seemingly 
lower amount of gas and dust in the inner envelope. If clumps 
are responsible for the gradient in the dust-to-gas ratio, we have 
to assume that the optical depth effect caused by clumping is 
more severe for dust than for gas. Considering that a cloud of 
gas particles experiences an internal pressure, whereas a cloud 
of dust particles does not, this could be a valid assumption. We 
note that a clumped wind is also invoked bv iDiikstra et al.l (l2006l) 
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to explain the observed high crystalline H2O ice fraction in OH 
127.8+0.0. 


5. Conclusions 

We have combined two state-of-the-art radiative transfer codes, 
MCMax for the continuum radiative transfer, and GAS- 
TRoNOoM for the line radiative transfer. We justified the use of 
more consistent dust properties in the gas modeling by showing 
that the dust component of the CSE has a significant infiuence 
on the excitation of H2O at high mass-loss rates, while the dust 
condensation radius is important for both CO and H2O at low 
mass-loss rates. 

We presented new PACS data of OH 127.8+0.0, the first 
AGB OH/IR star for which a far-IR spectrum was taken with 
this instrument. We applied our approach to the combination 
of the PACS spectrum, HIFI observations of two CO transitions 
taken in the framework of the SUCCESS Herschel Guaranteed 
Time Program, ground-based JCMT observations of low-J CO 
transitions, and the ISO-SWS and ISO-LWS spectra. The com¬ 
bination of the HIFI and ground-based observations suggests a 
discrepancy between the lowest-J (7=1-0 and 7 = 2-1) and 
the higher-J (7 = 3-2 and up) CO lines, which may point to a 
recent onset of a superwind in OH 127.8+0.0. The IR continuum 
is modeled with a dust composition of metallic iron, amorphous 
silicates, crystalline silicates (forsterite and enstatite), and amor¬ 
phous and crystalline H2O ice. We found a dust-mass-loss rate 
of Md = (5 + 1) X 10“^ Mo yr“^ and a contribution of H2O 
ice to the total amount of dust beyond the H2O ice condensa¬ 
tion radius of (16 + 2)% with a crystalline-to-amorphous ratio 
of 0.8 + 0.2. The CO transitions are modeled with an empir¬ 
ical temperature law resulting in four models with a constant 
gas-mass-loss rate ranging between Mg = 1.0 X 10-4 yj.-l 
and Mg = 0.2 X lO'^ Mq yr"*, accurate to within a factor of 
three. The older mass-loss episode, traced by the outer regions 
of the CSE, is estimated to be Mg ~ lx 10“^ M© yr“^ with 
the transition between the low and high mass-loss rate occurring 
at Rvm ~ 2500 - 4000 We derived a critical H2O vapor 
abundance of (1.7 + 0.2) x 10“^ from the H2O ice content of 
the CSE. This constrains the minimum amount of H2O vapor 
required to produce the observed amount of H2O ice assuming 
100% freeze-out efficiency. We note that the comparison be¬ 
tween H2O vapor models and the PACS spectrum shows a flux 
overestimation at shorter wavelengths and a flux underestima¬ 
tion at longer wavelengths. Even though these differences are 
within the absolute flux calibration of the PACS instrument, the 
wavelength-dependent discrepancy cannot be explained. 

We derived the dust-to-gas ratio following three methods, 
which are sensitive to different regions of the outflow. We found 
for the first time indications of a gradient in the dust-to-gas ratio 
with radial distance from the star. Possible explanations for this 
behavior can include dumpiness, variable mass loss, or contin¬ 
ued dust growth beyond the condensation radius, of which the 
first suggestion seems the most likely. 

Additionally, we reported the first detection in an AGB cir- 
cumstellar environment of OH cascade rotational lines involved 
in the far-infrared pumping mechanism of the 1612 MHz OH 
maser. 
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Appendix A: Integrated line strengths 

Table lA.ll lists integrated line strengths of all detected ortho- and 
para-H20 vapor emission lines and the 1612 MHz OH maser 
formation rotational cascade lines in the PACS spectrum shown 
in Figs. [T0| and [TT] Because the OH emission lines occur in 
doublets, the integrated lin e strengths for both com ponents have 
been summed. We refer to I Sylvester et al.l (1 19971) for details on 
OH spectroscopy. Where confusion due to line blending occurs, 
we indicate this clearly, as well as list all H2O transitions that 
may contribute to the emission line. As such we cannot distin¬ 
guish the relative contribution of each transition in the blend. 
Blends that might be caused by the emission of other molecules 
not modeled in this study are not indicated. 


Article number, page 17 of [19] 






Table A.l. Integrated line strength Fml (W/m^) of detected ortho- and para-H20 vapor emission lines and the 1612 MHz 
OH maser formation rotational cascade lines in the PACS spectrum shown in Figs. [T0| and [TT] The rest wavelength To 
(jum) of the transition is indicated. The percentages between brackets indicate the uncertainty on FiM, which includes 
both the fitting uncertainty and the PACS absolute fiux calibration uncertainty of 20%. 


PACS 

Molecule 

Vibrational 

Rotational 

Ao 

F int 

band 


state 

transition 

jum 

(W/m2) 

B2A 

0 -H 2 O 

y = 0 

‘^Ka,Kc - ^ 3,2 - 32,1 

58.70 

1.99e-16(25.8%) 


0 -H 2 O 

V3 = 1 

‘^Ka,Kc - ^ 3,1 - 32,2 

58.89 

6.64e-17(45.2%) 


P-H 2 O 

y = 0 

^Ka,Kc = ^2,6 - 61,5 

59.99 

* 1.51e-16(38.6%) 


P-H 2 O 

y = 0 

‘^Ka,Kc - ^2,6 - ^ 3,5 

60.16 

>’8.906-17(34.8%) 


0 -H 2 O 

y 2 = 1 

^Ka,Kc = 33,0 - 22,1 

60.49 

>’8.656-17(40.6%) 


0 -H 2 O 

y = 0 

^Ka,Kc - ^6,1 - 65,2 

63.91 

*4.506-16(23.1%) 


0 -H 2 O 

y = 0 

‘^Ka,Kc - ^6,0 - 65,1 

63.93 



0 -H 2 O 

y 2 = 1 

^Ka,Kc = 80,8 - 7 lJ 

63.95 



0 -H 2 O 

y = 0 

^Ka,Kc = 76,1 - 15,2 

63.96 



0 -H 2 O 

y = 0 

7 Ka,Kc = 62,5 - 5 1,4 

65.17 

*’>’ 1.316-16(35.2%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 33,0 - 22,1 

66.44 

>’ 1.286-16(29.0%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 33,1 - 22,0 

67.09 

*2.596-16(23.1%) 


P-H 2 O 

y 2 = 1 

7 Ka,Kc = ^ 2,4 - 4 i ,3 

67.26 

* 1.526-16(25.2%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 33,0 - 3 o ,3 

67.27 



0 -H 2 O 

y 2 = 1 

7 Ka,Kc = 32,1 - 2i,2 

70.29 

7.596-17(27.8%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 82,7 - 81,8 

70.70 

1.316-16(22.4%) 


P-H 2 O 

y = 0 

7 Ka,Kc = ^ 2,4 - 4 i ,3 

71.07 

>’ 1.926-16(22.6%) 


P-H 2 O 

y = 0 

7 Ka,Kc = ^ 1,7 - 60,6 

71.54 

1 . 116 - 16 ( 29 . 2 %) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 7 o ,7 - 61,6 

71.95 

1.286-16(26.3%) 

B2B 

0 -H 2 O 

y 2 = 1 

7 Ka,Kc - ^2,1 - 2 i,2 

70.29 

* 1.146-16(34.4%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 82,7 - 81,8 

70.70 

1.276-16(29.1%) 


P-H 2 O 

y = 0 

7 Ka,Kc - ^ 2,4 - 4 i ,3 

71.07 

>’ 1.886-16(22.8%) 


P-H 2 O 

y = 0 

7 Ka,Kc -^ 1,7 ~ 60,6 

71.54 

9.546-17(31.8%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 7 o ,7 - 61,6 

71.95 

1.476-16(25.8%) 


P-H 2 O 

y = 0 

7 Ka,Kc - ^ 3,7 - 92,8 

73.61 

9.096-17 (26.0%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 72,5 - 63,4 

74.95 

1.266-16(22.7%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 32,1 - 2i,2 

75.38 

>’ 1.056-16(23.7%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 85,4 - 84,5 

75.50 

5.406-17(31.8%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 55,1 - 54,2 

75.78 

*3.116-16(21.1%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 75,3 - 74,4 

75.81 



P-H 2 O 

y = 0 

7 Ka,Kc = 65,2 - 64,3 

75.83 



0 -H 2 O 

y = 0 

7 Ka,Kc = 55,0 - 54,1 

75.91 

* 1.096-16(26.8%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 65,1 - 64,2 

76.42 

1.086-16(25.4%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 75,2 - 74,3 

77.76 

>’7.396-17(25.6%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 42,3 - 3 1,2 

78.74 

1.436-16(22.8%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 61,5 - 52,4 

78.93 

*8.616-17(26.6%) 


P-H 2 O 

y 2 = 1 

7 Ka,Kc = 64,3 - 63,4 

78.95 



OH 

y = 0 

2 ni/ 2 / = i /2 -2 n3/2/ = 3/2 

79.12-79.18 

‘’3.096-16 (35.9%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 83,6 - 82,7 

82.98 

9.186-17(25.7%) 


P-H 2 O 

y 2 = 1 

7 Ka,Kc = 32,2 “ 2i,i 

83.24 

5.266-17 (39.8%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 60,6 - 51,5 

83.28 

8.806-17(28.1%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 7 i ,6 - 7 o ,7 

84.77 

* 1.566-16(22.9%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 84,5 - 83,6 

85.77 

* 6.926-17 (25.5%) 


0 -H 2 O 

y 2 = 1 

7 Ka,Kc = 64,2 - 63,3 

85.78 



P-H 2 O 

y = 0 

7 Ka,Kc = 32,2 “ 2i,i 

89.99 

* 1.496-16(22.9%) 


P-H 2 O 

y = 0 

7 Ka,Kc = 74,4 - 73,5 

90.05 

8.096-17 (24.0%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 64,3 - 63,4 

92.81 

>’ 1 . 206 - 16 ( 21 . 8 %) 


P-H 2 O 

y = 0 

7 Ka,Kc = 73,5 - 72,6 

93.38 

9.906-17 (22.4%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 65,2 - 72,5 

94.17 

* 1.076-16(23.7%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 54,2 - 53,3 

94.21 



0 -H 2 O 

V 3 = 1 

7 Ka,Kc - 74,4 - 73,5 

94.61 

* 1.106-16(23.4%) 


0 -H 2 O 

y = 0 

7 Ka,Kc = 62,5 - 61,6 

94.64 



0 -H 2 O 

y = 0 

7 Ka,Kc - 44,1 - 43,2 

94.71 

1.406-16(22.5%) 


Notes. Line strengths fiagged for potential line blends (see Sect. 12.2.11 . Transitions that might cause the line blend 
are mentioned immediately below the flagged transition. Transition detected in an emission doublet. The given value 
is the sum of the line strengths of both emission lines in the do ublet. The selection of H 2 O vapor emission lines based 
on which we have derived the H 2 O vapor abundance (see Sect. 14.2.31 ). 
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Table A.l. continued. 


PACS 

Molecule 

Vibrational 

Rotational 

Ao 

^ int 

band 


state 

transition 

jum 

(W/m^) 


P-H 2 O 

V 2 = 1 

^Ka,Kc = 51,5 - 4 o,4 

94.90 

4.91e-17(32.5%) 


0 -H 2 O 

y = 0 

^Ka,Kc = 94,5 - 85,4 

95.18 

2.26e-17 (40.0%) 


P-H 2 O 

y = 0 

‘^Ka,Kc = 51,5 - 4 o,4 

95.63 

8.59e-17(28.8%) 


P-H 2 O 

y = 0 

^Ka,Kc = 44,0 - 43,1 

95.88 

'’9.46e-17(26.3%) 


0 -H 2 O 

y = 0 

‘^Ka,Kc = 54,1 - 53,2 

98.49 

*6.85e-17(51.4%) 


OH 

y = 0 

^Uif2J = 5/2 - 2 ni/ 2 / = 3/2 

98.72-98.74 

*’M.54e-16(29.6%) 


P-H 2 O 

V 2 = 1 

‘^Ka,Kc = 83,5 - 74,4 

98.73 


RIA 

P-H 2 O 

y = 0 

7Ka,Kc = 64,2 - 63,3 

103.92 

^1.70e-16(23.4%) 


P-H 2 O 

y = 0 

7Ka,Kc = 61,5 - 60,6 

103.94 



0 -H 2 O 

y = 0 

7Ka,Kc - ^3,4 - 62,5 

104.09 

1.27e-16(26.4%) 


0 -H 2 O 

V 3 = 1 

7Ka,Kc - ^ 2,0 - ll,l 

104.81 

4.15e-17(39.1%) 


0 -H 2 O 

y = 0 

7Ka,Kc = ^ 2,1 - ll ,0 

108.07 

M.25e-16(21.4%) 


P-H 2 O 

y = 0 

7Ka,Kc - ^2,4 - 5 1,5 

111.63 

9.51e-17(24.1%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 74,3 - 73,4 

112.51 

3.53e-17 (33.9%) 


0 -H 2 O 

V 2 = 1 

7Ka,Kc = 64,3 - 7 1,6 

112.80 

*3.20e-17(39.4%) 


0 -H 2 O 

y = 0 

7Ka,Kc - 44,1 - 51,4 

112.80 



0 -H 2 O 

V 3 = 1 

7Ka,Kc = ^2,4 - 5 1,5 

112.89 



0 -H 2 O 

y = 0 

7Ka,Kc = 4 i,4 - 3o,3 

113.54 

9.06e-17 (23.7%) 


P-H 2 O 

y = 0 

7Ka,Kc = 53,3 - 52,4 

113.95 

1.24e-16(21.9%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 43,2 - 42,3 

121.72 

^9.60e-17(21.2%) 


P-H 20 

y = 0 

7Ka,Kc = 84,4 - 83,5 

122.52 

2.63e-17 (34.0%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 93,6 - 92,7 

123.46 

3.11e-17(30.1%) 


0 -H 2 O 

V 2 = 1 

7Ka,Kc = 51,4 - 5 o,5 

124.85 

^5.90e-17(24.5%) 


P-H 2 O 

y = 0 

7Ka,Kc = 4 o,4 - 31,3 

125.35 

7.52e-17 (22.4%) 


P-H 2 O 

y = 0 

7Ka,Kc = ^3,1 - 32,2 

126.71 

^8.530-17(21.6%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 72,5 “ 7 i,6 

127.88 

*6.84e-17(26.4%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 94,5 - 93,6 

129.34 

4.27e-17 (24.0%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 42,3 - 4 i,4 

132.41 

8.39e-17(21.3%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 51,4 - 5 o,5 

134.94 

^7.45e-17(22.2%) 


0 -H 2 O 

y = 0 

7Ka,Kc = ^3,0 - 32,1 

136.50 

8.01e-17(21.7%) 


P-H 2 O 

y = 0 

7Ka,Kc = 31,3 - 2 o,2 

138.53 

*7.59e-17(21.7%) 


P-H 2 O 

y = 0 

7Ka,Kc = 84,4 - 75,3 

138.64 



P-H 2 O 

V 2 = 1 

7Ka,Kc = 63,3 - 62,4 

140.06 

*2.99e-17(31.7%) 


P-H 2 O 

y = 0 

7Ka,Kc = 4 i,3 - 32,2 

144.52 

M.63e-17 (23.5%) 

RIB 

P-H 2 O 

y = 0 

7Ka,Kc - 4 i,3 - 32,2 

144.52 

''4.806-17(23.7%) 


P-H 2 O 

y = 0 

7Ka,Kc = 43,1 - 42,2 

146.92 

6.726-17(22.0%) 


0 -H 2 O 

V 3 = 1 

7Ka,Kc - ^ 3,5 - 82,6 

148.64 

*5.126-17(25.3%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 83,5 - 82,6 

148.71 



0 -H 2 O 

y = 0 

7Ka,Kc - ^4,2 - 61,5 

148.79 



0 -H 2 O 

V 2 = 1 

7Ka,Kc - ^ 2,1 - 2 i ,2 

153.27 

1.516-17(33.8%) 


P-H 2 O 

y 2 = 1 

7Ka,Kc - ^2,4 - 61,5 

154.02 

1.016-17(45.7%) 


P-H 2 O 

y = 0 

7Ka,Kc = ^ 2,2 - 3 1,3 

156.19 

* 1.286-16(20.7%) 


P-H 2 O 

y = 0 

7Ka,Kc - ^2,3 - 43,2 

156.27 



0 -H 2 O 

y = 0 

7Ka,Kc = 53,2 - 52,3 

160.51 

4.606-17 (23.6%) 


OH 

y = 0 

2ni/2/ = 3/2 -2ni/2/= 1/2 

163.12-163.4 

1.476-16 (30.5%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 73,4 - 72,5 

166.81 

*3.776-17(31.7%) 


0 -H 2 O 

V 3 = 1 

7 Ka,Kc = 62,4 - 61,5 

166.83 



P-H 2 O 

y = 0 

7 Ka,Kc = 62,4 - 61,5 

167.03 

4.116-17(27.9%) 


P-H 2 O 

y = 0 

7Ka,Kc = 63,3 - 62,4 

170.14 

*4.506-17(35.1%) 


P-H 2 O 

y = 0 

7Ka,Kc = 53,3 - 60,6 

174.61 

*8.256-17(21.9%) 


P-H 2 O 

y = 0 

7Ka,Kc = 3 o,3 - 2i,2 

174.63 



P-H 2 O 

V 3 = 1 

7Ka,Kc = 3 o,3 - 2i,2 

174.66 



0 -H 2 O 

y = 0 

7Ka,Kc = 43,2 - 5o,5 

174.92 

1.496-17(43.1%) 


0 -H 2 O 

y = 0 

7Ka,Kc = 2 i ,2 - lo,l 

179.53 

4.936-17 (24.1%) 


0 -H 2 O 

y = 0 

7Ka,Kc = ^2,1 - 2i,2 

180.49 

''5.876-17(26.7%) 


Notes. Line strengths flagged for potential line blends (see Sect. 12.2.11 . Transitions that might cause the line blend 
are mentioned immediately below the flagged transition. Transition detected in an emission doublet. The given value 
is the sum of the line strengths of both emission lines in the doublet. The selection of H 2 O vapor emission lines based 
on which we have derived the H 2 O vapor abundance (see Sect. l4.2.3l l. 
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